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1. INTRODUCTION 



1.1 PURPOSE . 

The purpose of this investigation is to show how adjustment compu- 
tations may be applied to certain coastal hydrographic surveying opera- 
tions. Particular emphasis is placed upon the determination of shore 
positions using information obtained aboard the survey vessel from vis- 
ual and electronic positioning systems. The inverse problem, determi- 
nation of the position of the ship by means of observations made at 
known shore positions, is also treated with an adjustment computation 
procedure. 

The presently employed visual and electronic methods of location 
in coastal hydrographic surveying are described in order to illustrate 
the particular problems to which the adjustment procedures are to be 
applied. Tlie instruments and equipment currently utilized to obtain 
positioning information are explained to indicate their relative merits 
and the accuracies that can be expected with them. 

A brief summary of the state of automation in hydrographic sur- 
veying is presented to illustrate the feasibility of incorporating 
automated adjustment procedures into these operations on a real-time 
basis. 

After this background information has been presented, the adjust- 
ment procedures are described for the two general problems. In each 
instance the overall adjustment procedure to be employed is outlined 

using standard matrix notations. The use of matrix algebra was decided 
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upon because of its straight forward adaptation to high-speed digital 
computer operations. Next, the various equations and relationships 
pertaining to the particular problem being discussed are developed in 
detail. The error analysis method for each problem is then explained. 
Finally, a numerical example of each adjustment is presented to illus- 
trate the application of the procedure to a hydrographic survey opera- 
tion. These numerical examples are solved by computer programs which 
are designed to determine all necessary auxiliary information and per- 
form the required adjustment and error analysis computations. Partic- 
ular attention has been given to developing the programs in as general 
a form as possible to allow for a wide latitude in the type and amount 
of observational input information. 

Since real input data could not be obtained, simulated information 
was compiled for each problem using plotting sheets and equipment pro- 
vided by the U.S. Naval Oceanographic Office. For this reason the ac- 
curacies of determinations must not be construed as being representa- 
tive of results which would be obtained with the use of actual observa- 
tional data. The primary objective of the numerical examples is to il- 
lustrate the methods in which the adjustment computations would be ap- 
plied to coastal survey problems. 

1.2 APPLICATIONS . 

1.2.1 Shipboard Observation Application . 

A procedure is developed which permits the determination of ship 
and shore positions using only shipboard information. The data pro- 
vided by the solution of this problem are utilized for charts to aid 
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in the conduct of navigation and offensive military operations along 
inaccessible coastal regions. This application is designed in the form 
of a rapid self-contained positioning system which will enable the sur- 
vey vessel to locate itself with relation to the coastal topography 
without placing personnel or equipment ashore. It could find particu- 
lar use as a method for conducting surveys in hostile environments. As 
such, this procedure represents a possible method for arriving at a so- 
lution to certain special coastal hydrographic survey problems de- 
scribed by THOMAS (1). 

The problem is dealt with in a completely generalized form in this 
investigation. That is, it is assumed that no topographic survey data 
is available concerning the positions of terrestrial objects. Further, 
it is assumed that the coast is inaccessible, either by virtue of im- 
penetrable natural obstacles or, more likely, because the territory is 
occupied by an unfriendly nation. A second objective of the hydro- 
graphic survey then, in addition to obtaining the usual ocean depth 
information, is to produce data for use in developing a rough topo- 
graphical survey; to the extent of determining the positions of promi- 
nent natural and man-made objects. These terrestrial points would la- 
ter be used by combatant ships for coastal navigation and as references 
for offshore positioning preparatory to offensive bombardment opera- 
tions against coastal targets. 

As a modification to the completely generalized solution, the sit- 
uation frequently arises wherein limited positional information is 
available concerning terrestrial objects. This information, combined 
with the shipboard observations, could be utilized to determine the 
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locations of newly constructed buildings such as enemy fortifications, 
as well as the usual ship positions, with considerably more accuracy 
than would be possible with the completely generalized problem. 

1.2.2 Shore-based Observation Application . 

The primary application of the ship positioning information pro- 
vided by visual shore observations is for the precise location of 
soundings. These data are subsequently utilized for the compilation 
of accurate depth information on nautical charts. Other uses of this 
method of accurate ship positioning include geophysical surveys, plan- 
ning for offshore military defense and civil engineering projects, lo- 
cation of submerged objects, cable laying operations, harbor improve- 
ment planning and oceanographic studies. 

Since this method yields redundant information, an adjustment is 
performed to obtain the best possible ship position. The results ob- 
tained for ship positions along a considerable length of the coast, 
together with the positions as determined by the ship itself simulta- 
neously from an electronic system, can provide the information required 
to calibrate and check the accuracy of the electronic system. Such 
comparisons could disclose systematic errors in electronic lattices 
along the coast which would reduce their accuracy. 0SB0RN(2) discusses 
in detail this method of determining the accuracy of and calibrating 
an electronic surveying system at the time of initial installation. 
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2. HYDROGRAPHIC SURVEYING 



2.1 GENERAL . 

Hydrographic surveying consists of the measurement of ocean depths 
and determination of the nature of the sea floor. The locations of 
certain features of adjoining coastal areas are also ascertained for 
use in aiding coastal navigation. The hydrographic survey can include 
gravimetric, tidal, current and meteorological data measurements as 
well as the usual position and sounding information. 

The positions at which depth measurements are made are usually de- 
termined by reference to established points on shore. This procedure 
enables the offshore points to be located at their correct geographic 
positions, and insures that the land and marine features are in proper 
relationship to each other. The principal object of obtaining water 
depth and bottom composition and configuration information by hydro- 
graphic surveying is to enable the compilation of nautical charts and 
related publications for the mariner(3) . 

Depth or sounding information is obtained almost exclusively with 
the electronic echo sounder. Only in extremely shallow and inaccessi- 
ble water areas are lead lines and sounding poles used in place of echo 
sounding devices. Echo sounding equipment of the recording type is 
usually used since it provides a graphic profile of the ocean floor, 
enabling the detailed charting of submarine relief. 

2.2 COASTAL HYDROGRAPHY . 

This investigation is concerned with hydrographic survey opera- 
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tions which are accomplished along the coast in sight of land. This 
branch of hydrographic surveying involves the determination of the 
shoreline, limited coastal topographic features, and the ocean bottom 
configuration in the immediate vicinity of the coast. The marine re- 
gion encompasses areas extending to ranges of ten to fifteen miles from 
the shoreline. 

Normally, prior to the conduct of the coastal hydrographic survey, 
a geodetic control survey is accomplished in order to establish control 
for shore positions. These shore points are used either to position 
target signals for visual sighting from the survey ship or for the lo- 
cation of sites for electronic positioning system transmitters. The 
usual interval between signals for coastal surveys is one to two miles. 
The distances may be reduced to one-half to one mile for large scale 
harbor surveys. 

Small launches, known as sounding boats, are utilized to obtain 
depth measurements close to shore. In areas where the bottom slopes 
gradually, systems of sounding lines are normally taken parallel to the 
coast. The spacing between lines is least nearest the shoreline, usu- 
ally less than fifty yards. It is increased outward from the shore. 

The spacing depends primarily on the depth of water and the type of 
bottom. 

In addition to the survey of coastal waters, information is ob- 
tained on the shoreline itself and immediately adjacent land areas and 
outlying islands. This coastal mapping is generally limited to terrain 
features needed for control of the survey. Conspicuous landmarks which 
can be used by the mariner as aids in alongshore navigation are also 
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accurately located for inclusion on the nautical chart(4) . 

The most economical method for obtaining the preliminary shoreline 
and coastal topographic data, when it is feasible, is by means of aeri- 
al photography. Large-scale photogrammetric shoreline manuscript maps 
are produced to furnish advance reconnaissance information and provide 
detail to aid the hydrographer in mapping the coastline and terrestrial 
features. 

2.3 VISUAL POSITIONING . 

Various procedures are utilized to accurately determine the posi- 
tion of the survey craft from which depth and other measurements are 
made. The two general methods of hydrographic control most commonly 
used are visual and electronic. 

Visual control is normally the method employed when surveying is 
performed within sight of the coast. Visual reference is made to con- 
spicuous landmarks; objects such as geodetic survey control signals, 
prominent natural or man-made objects and supplementary stations estab- 
lished especially for the hydrographic survey(see Figure 2.1). The 
supplementary signals are positioned in relation to geodetic control 
stations by intersection, resection or traverse methods. 

The positioning of a survey vessel by visual means usually falls 
within one of the three following categories: 

(1) Resection from the ship. 

(2) Intersection from the shore. 

(3) Simultaneous intersection-resection. 

A number of procedures are available to the hydrographic surveyor by 
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Figure 2.1, Hydrographic signal erected at a known shore 
point for use in controlling a coastal survey. 
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which he may apply the three basic methods. The requirements of the 
particular survey, primarily the resources available for its conduct 
and the accuracy specified for the results, dictate the particular pro- 
cedure to be employed. The environment of the area to be surveyed also 
influences the specific methods to be utilized in accomplishing the 
project. 

Both of the problems presented in this investigation depend upon 
visual observations for an integral part of their input information. 
Also, both utilize certain of the basic concepts of visual location 
methods for their solution. For these reasons and in order to provide 
suitable background information for the investigation, certain visual 
location and solution methods are described below. 

2.3.1 Location Methods . 

(a) By two angles from the ship. 

This is the most common method used in conventional coastal sur- 
veying operations. Visual reference is made to conspicuous landmarks 
of known location along the shore. Usually two horizontal sextant an- 
gles, designated a and 6 in Figure 2.2, are observed simultaneously 
from the ship in the manner indicated in Figure 2.3. If conditions 
permit, a third "check angle", the angle 0 subtended by the extreme 
points T and V in Figure 2.2 ( is also observed to assist in detecting 
gross errors. A less accurate observation method entails the use of 
compass bearings to the shore objects instead of sextant angle cuts. 

In practice, a mechanical solution for the resultant ship's position 
is obtained with a device known as the three-arm protractor. This pro- 
cedure as well as the graphical and analytical solutions are described 
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Figure 2.2. Location by two angles at the survey ship. 



in Article 2.3.2. 

The three-point resection locating method is advantageous from the 
standpoint that the entire operation takes place aboard the survey 
craft. Also, if a sufficient number of well-defined prominent objects 
of known location are visible along the coast, the task of establishing 
preliminary shore control is considerably reduced. A significant dis- 
advantage, and it applies to all conventional visual positioning pro- 
cedures, lies in the fact that they may be performed only in daylight 
hours and during periods of good visibility. 

Another serious disadvantage of the resection method results from 
the fact that the ship is a continuously moving observation platform. 

In addition, the instruments are of low accuracy or have relatively 

slow readout capabilities. Hence, when greatest accuracy is required 
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Figure 2.3. Observing sextant angles to prominent points ashore from 



the intersection method must be used, 

(b) By directions from shore. 

This method consists of determining the position of the survey 
craft by simultaneously observing intersecting directions to the ship 
from two or more control points ashore. The shore stations are nor- 
mally determined in advance by triangulation, traverse or trilatera- 
tion. 

With the ship either moored or drifting slowly, the simultaneous 
directions are observed to it from the shore stations, denoted as T, 




Figure 2.4. Location by directions from shore stations. 

U, V and W in Figure 2.4, with either a theodolite or a continuously 
tracking azimuth instrument. A radio system is employed to coordinate 
the observation procedures and to relay direction data to the ship for 
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plotting. 



The principal advantage of the intersection method is that obser- 
vations are made from accurately determined points. For this reason 
and because the observing platform is stable, higher precision instru- 
ments may be employed. These factors result in a much higher accuracy 
in the ship position location than is possible with the resection meth- 
od from the vessel. Another advantage of this method is that the pre- 
liminary task of setting out and erecting shore signals is avoided. 
Also, the problem of maintaining the vessel in a precise location in 
areas where currents are strong is eliminated. Disadvantages accrue 
from the separation of the survey party members and the associated 
problems which can occur with faulty radio communications between the 
ship and shore stations. 

(c) By ship angles and shore directions. 

This procedure, the simultaneous resection-intersection problem, 
is essentially a combination of methods (a) and (b) . Since the obser- 
vations taken from the shore stations are considerably more accurate 
than the ship information, this method serves primarily only as a check 
on the ship observations. In practice it is seldom used, 

2.3.2 Resection Solution Procedures . 

Three basic methods are available for determining the sextant fix 
position. They are all examples of the solution of the classic three- 
point resection problem. Using the notations in Figure 2.2, the prob- 
lem may be stated as follows: given three shore points T, U and V and 
the values a and 6 of the observed sextant angles at the ship; to de- 
termine the position S of the ship. As mentioned in Article 2.3.1, the 
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Figure 2.5. Plotting sextant angle information aboard a survey vessel with a three- 
arm protractor for mechanical solution of the resection problem. 



procedure most widely used in practice is the mechanical solution shown 
in Figure 2,5. Two others, the graphical and the analytical, are also 
available. However, for practical reasons, primarily the excessive 
time required, they are not usually employed in hydrographic survey 
operations. 

(a) Mechanical solution. 

In this procedure, the two horizontal sextant angles a and 6 are 
set into the three-arm protractor. The three arms of the protractor 
are next made to coincide with the three shore objects T, U and V which 
have been transferred to the plotting sheet. The center of the pro- 
tractor then indicates the position S of the survey craft at the time 
of the observation, A small hole located at the center of the protrac- 
tor permits the position of the fix to be marked on the plotting sheet. 

(b) Graphical solution. 

Several methods are available by which a graphical solution to the 
resection problem may be obtained. Three such constructions are de- 
scribed in (5). A representative procedure is outlined below. Refer- 
ring to Figure 2.6, the lines between T and U and between U and V are 
first joined. At T and U the lines TOj and UO^ respectively are set 
off from TU as indicated with angles equal to 90° minus the observed 
sextant angle a which was subtended at S between T and U. From the 
line UV in a similar manner, the lines UC^ and VO^ are set off at 90° 
less the sextant angle 6. Two circles are then drawn, the first with 
center at 0^ and passing through T and U, and the second with center at 
0-, and passing through U and V. The desired position of the ship is 
the intersection S of the two circles. That this is so may be shown by 
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Figure 2.6. Graphical solution of the three-point problem, 
the fact that the angle at S between T and U at the circumference of 
circle 0^ is equal to half the angle at Oj between T and U on the same 
chord. But this angle is 2a. Thus the angle at S between T and U is a. 
In similar fashion it can be shown that the angle at S between U and V 
is 6. 

(c) Analytical solution. 

Numerous procedures have been developed for the analytical solu- 
tion of the three-point problem. A discussion of two such methods and 
an extensive theoretical and practical treatment of sextant positioning 
is given in (6). The following procedure is taken from (5). Referring 
to Figure 2.7, the lengths TU and UV and the angles a, 6 and y are con- 
sidered as the given quantities. Using the law of sines, 

SU * . Py.,gf». f (2.1) 

sin a sin 8 

Now letting d ■ e + f (2.2) 



16 





equation (2,1) can be written as: 

SU « ^ sin 6 = W sin(d-e) 
sin a sin 6 

Using the relationship for the sine of the difference of two angles, 

the right side of (2,3) becomes: 

TU sin e UV(sin d cos e - cos d sin e) 

sin a STTb L (2 * 4) 

Solving (2,4) for e and simplifying, 

e * ^‘ [nHlnV sin J * d ] < 2 - 5 > 

With e from (2,5), the angle f can be computed from (2,2) since d is a 

known quantity, given by: 

d * 360° - (o ♦ 8 ♦ y) (2.6) 

Finally, the distances SU, ST and SV are calculated using the law of 

sines. 
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2.3.3 Intersection Solution Procedures. 



As in the case of the resection problem, numerous acceptable me- 
thods are available for the solution of the intersection problem. The 
graphical and mechanical procedures consist simply of extending direc- 
tions or plotting angles from fixed or known positions. The intersec- 
tions formed by the known direction lines determine the location of the 
desired points. 

An analytical solution of the intersection problem is described in 
Chapter 5. The problem is explained in detail in conjunction with the 
associated adjustment computation. For this reason it is not treated 
separately in this chapter. 

2.3,4 Equipment Characteristics . 

(a) Hydrographic sextant. 

The primary instrument employed in the visual methods of hydro- 
graphic surveying is the hand-held hydrographic sextant which is used 
to measure the horizontal angle between two shore objects. It has a 
maximum range of about 140° and is capable of measuring an angle to the 
nearest minute. The principle of the sextant is that the deviation of 
a ray of light reflected successively from two mirrors is twice the an- 
gle between them. One mirror, a large index glass, is mounted on an 
axis about which it can be rotated. The variable angle between the in- 
dex glass and the fixer 1 horizon glass mirror formed when the two ob- 
jects sighted are brought into coincidence measures the angle subtended 
by the objects. 

There are several sources of error in the sextant, some of which 
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may be adjusted for by the user and some which may not. 0SB0RN(2) 
gives a complete description of the various errors and how they may be 
minimized or removed. The four principal non-adjustable errors are 
summarized below. 

(1) Prismatic error which occurs when the mirrors are not 

parallel. 

(2) Graduation errors in the arc, micrometer drum and ver- 
nier of the sextant which are present when the instrument has been im- 
properly constructed or calibrated. 

(3) Centering error which results when the index arm is not 
pivoted at the exact center of curvature of the arc. 

(4) Telescope error, resulting from resolving power, magni- 
fication and curvature of field characteristic limitations. 

The four adjustable errors which together comprise the index error 

are : 



(1) The error which results when the index mirror has not 
been adjusted so as to be perpendicular to the frame of the sextant. 

(2) Side error which occurs when the horizon glass is not 
adjusted so as to be perpendicular to the frame of the sextant. 

(3) Nonparallelism of the index mirror and the horizon glass 
when the index arm is set exactly to zero. 

(4) Nonparallelism of the telescope and the frame of the 

sextant . 

The determination of the errors of a sextant is normally based on 
a series of observations. With careful adjustments and proper hand- 
ling, the sextant will remain in aujustment for a considerable length 
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of time. Checks are made, however, prior to and after each operation. 

The greatest sources of error in sextant angle observations result 
from the shore objects being observed. This is particularly so when 
the shore points are buildings, water towers and natural objects which 
do not present a clearly defined narrow vertical line for sighting. 

Only with such objects as radio antennas and shore signals constructed 
by the survey party can maximum accuracy be attained. Poor contrast 
between the signal and the background region can also seriously degrade 
the accuracy of an observation. Another factor which normally far out- 
weighs instrument limitations is the continuous and usually irregular 
movement of the survey vessel. This rolling and pitching can severely 
limit the observer's ability to obtain an accurate sighting. Indeed, 
observations must often be discarded and repeated when an observer is 
unable to make a proper sighting at the designated fix time. 

Atmospheric conditions which cause haze, glare and similar pheno- 
mena; and the inability of the human eye to detect objects that subtend 
an angle smaller than a minute of arc also contribute to reduced accu- 
racies. 

It should be clear from the foregoing discussion that it is almost 
impossible to assign a meaningful standard accuracy value to a sextant 
angle observation. The conditions which exist during the survey must 
be analyzed and weighed in order to arrive at satisfactory values to be 
assigned for a particular project. 

f ) Ship's gyro compass. 

Compass bearings are much less preferable than sextant angle cuts 
in normal coastal hydrographic surveying operations as indicated in 
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Article 2.3.1. However, in the problem to be dealt with in Chapter 4, 
compass bearings must be taken to terrestrial objects for use as ob- 
served azimuths in the figural adjustment computations. For this rea- 
son, the error sources which influence and affect the resultant accu- 
racy of a gyro compass bearing are briefly described. 




Figure 2.8. Gyro compass direction components. 

Referring to Figure 2.8, C is the meridian convergence and E the 
gyro compass orientation error. The observed gyro compass bearing to 
the shore object is a^, while t^ is the observed plane azimuth of the 
shore object referenced to the grid system being employed. Thus t^ may 
be expressed as follows: 



t, ~ a + E + C (2.7) 

b g 

The gyro compass orientation error is the resultant of several 
systematic errors, certain of which are eliminated or compensated for 
in the construction of the compass and others which are minimized by 
adjustments made aboard ship. The five principal errors are enumerated 
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below. Detailed descriptions of these errors with procedures for elim- 
inating or minimizing each of them are given in B0WDITCH(3). 

(1) Speed Error. This error is introduced by the motion of 
the ship along its track. The direction of ship travel affects the 
magnitude of the error. It is corrected for in the design of the com- 
pass . 

(2) Damping Error. This error, sometimes known as the bal- 
listic damping error, is present in certain types of compasses. Its 
extent depends upon the method in which damping is accomplished in the 
compass. It is removed with suitable corrections. 

(3) Ballistic Deflection Error. This error occurs because 
an accelerating force acts on the compass causing a surge of mercury 
from one part of the system to another each time the north-south compo- 
nent of the ship's speed changes. It is minimized with certain compass 
adjustments. 

(4) Quadrantal Error. This error is caused by the swinging 
of the compass, an effect which varies according to the ship's heading. 
It is corrected for by adding weights to balance the compass so that 
the weight is the same in all directions from the center. 

(5) Gimballing Error. This error results when compass di- 
rections are measured in a tilted plane. When the compass is tilted, 
the projection of its outer rim onto the horizontal is an ellipse. 

Thus the graduations are not equally spaced with respect to a circle. 
The best way to avoid this error is to read the compass only when it 
is horizontal or nearly so. 

The residual gyro compass orientation error which remains after 
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all corrections and compensations have been applied, normally amounts 
to only about a tenth of a degree. It is applied equally to gyro bear- 
ings in all directions in order to obtain true bearings. 

Since the gyro compass can be read only to tenths of degrees, it 
can be assumed for the determination of azimuth values that the plane 
azimuth t and the projected geodetic azimuth T are the same. This is 
permissable because the (t-T) correction seldom exceeds several seconds 
of arc. 

A correction for the meridian convergence C must be computed and 
applied to convert true azimuths to grid azimuths if an X-Y coordinate 
system such as UTM is used for plotting. With a geographic latitude 
and longitude selected from the center of the project area plotting 
sheet, sufficient accuracy for computation of the difference between 
the true and grid azimuths is obtained with the following formula: 

d = AX sin $ (2.8) 

In this expression d is the difference in azimuths, AX is the differ- 
ence in longitude between the point and the central meridian and 4> is 
tne latitude of the point. The sign of the correction depends upon the 
hemisphere and the side of the central meridian on which the computa- 
tion point is located. 

The considerations in determining the accuracy of a gyro bearing 
direction observation are essentially the same as previously discussed 
in the case of the sextant angle observation. However, since the com- 
pass reading is limited to tenths of degrees and because the uncertain- 
ty in determining the residual compass orientation error may amount to 
almost a tenth of a degree, the resultant direction observation 



23 



accuracy limiting factor is caused by the instrument rather than by 
external influences. A series of observations is taken with the com- 
pass to ascertain the reading accuracy. This value is then combined 
with the instrument accuracy determined for the individual compass to 
arrive at an overall observation accuracy. 

2.3.5 Position Accuracy . 

Much has been written concerning the position accuracies which can 
be obtained with sextant angle measurements. A foremost requirement 
for obtaining satisfactory results is that the figure formed by the 
shore control points and the ship should be of maximum strength. Cer- 
tain general rules may be applied concerning the geometric strength of 
the resection figure. A fix will be weak in the case where the point 
to be located is near a circle which passes through the three known 




Figure 2.9. Geometric strength considerations in the resection 
problem. 
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points. Indeed, on this "danger circle" as it is known, the position 



(Sj) is indeterminate (see Figure 2,9), A fix will also lack strength 
if the new point (S^) is so far from the fixed points that the fixed 
base lengths between known points are ineffective. Conversely, if the 
resected point (S^) is near one of the known points it's position will 
be well determined. Any point (S^) inside the shaded triangle formed 
by the three fixed points will be strongly determined. 

With these and certain other less obvious criteria, the following 
rules have been formulated as the conditions for obtaining the best 
position when the resection principle is applied to the coastal hydro- 
graphic survey situation. The strongest fix occurs when, according to 
JEFFERS(7) : 

(1) The observer is inside the triangle formed by the three 
shore objects. 

(2) The three shore objects are in a straight line or the 
center object lies between the ship and the line joining the other two. 

(3) The sum of the two sextant angles is greater than 50° 
and preferably when neither angle is smaller than 30°, 

(4) Two of the shore objects a considerable distance apart 
are in range and the angle to the third is at least 45°. 

(5) At least one of the angles changes rapidly as the ship 
moves along it's track, 

(6) The distance between the center object and the left and 
rignt-hand objects is longer than the distance from the observer to the 
center object. 

A method devised for the solution of the three-point hydrographic 
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survey problem and the error of position that results from a known 
error in the sextant reading is described in detail by ALEXANDER(6) . 



2.4 ELECTRONIC POSITIONING . 

Electronic control is the primary location determining method em- 
ployed beyond the normal visibility from shore. Electronic control is 
also employed in coastal surveying during periods of low visibility and 
when adequate visual control is not available. Numerous electronic 
positioning systems are currently in use for hydrographic surveying 
control. Practically all of them are based upon either the circular 
(direct ranging) or hyperbolic method of line positioning. A few are 
composite systems, utilizing certain aspects of both basic methods. 

The principles and procedures for using the various methods are beyond 
the scope of this investigation. Numerous publications treat the vari- 
ous system in detail; among them LAURILA(8) and S.P. 39 of the Inter- 
national Hydrographic Bureau(9). For the purposes of this investiga- 
tion only a brief summary of the basic characteristics and relative 
merits of the systems need be mentioned. 

2,4.1 Circular Systems . 

The circular or direct ranging methods depend for positioning de- 
termination on the measurement of the round trip travel time of radio 
waves from the ship’s transmitter-receiver to accurately positioned 
shore stations. Families of concentric circles are produced as shown 
in Figure 2.10. The intersection of any two lines of position deter- 
mines the position of the vessel. 
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Since such systems require radio transmission by the ship, the 
number of users is limited. Time-sharing procedures must be employed 
when more than a single ship is to use the same shore equipment and 
radio frequencies at any particular time. Ranges of such systems suit- 
able for hydrographic surveying operations are generally limited to 
less than one hundred miles, due primarily to the high frequencies 
used. The all-weather capability of this and other electronic posi- 
tioning systems renders them much more suitable for most surveying work 
than the visual methods. 

The network geometry of the circular systems permits increased 
accuracy in positioning over larger areas in the grid compared with 
the hyperbolic systems because the spacing between curves of the same 
family remains constant. This characteristic is discussed in Article 
2 . 4 . 4 . 
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2.4,2 Hyperbolic Systems. 



The hyperbolic systems depend upon the measurement of the time or 
phase difference of radio waves emitted from synchronized transmitters 
located at known positions ashore. The shipboard receiver determines 
the difference in arrival of the radio signal pulses. These difference 
measurements for each set of two transmitters cause families of hyper- 
bolic lines of position to be generated with the transmitters at the 
foci. A minimum of two such systems are normally established as a 
unit, with a master transmitter acting as the control for two end or 
slave stations. The master station triggers the slave stations' pulses 
from it's own transmission. A typical hyperbolic network pattern is 
shown in Figure 2.11, 

As with the circular systems, the intersection of two or more 




Figure 2,11. Hyperbolic electronic positioning system 
network . 



28 



lines of position determines the location of the ship. The accuracy 
of positioning varies considerably within the network, depending on the 
angle of intersection formed between the hyperbolic curves and the sep- 
aration between curves of the same family at the location of the ship. 
Position accuracies and their determination are described in detail in 
Article 2.4.4, 

Since the survey craft requires only a receiver, the number of 
users of hyperbolic systems is unlimited. The lower frequencies uti- 
lized in these systems permit greater ranges than with the circular 
methods. 

2.4.3 Composite Systems . 

Ranging and hyperbolic systems may be combined into composite 
positioning systems. Such systems are normally developed to improve 
the geometry of the electronic lattice. They are usually constructed 
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to satisfy a particular requirement, often when the area to be surveyed 
does not lend itself to the satisfactory installation of either a rang- 
ing or a hyperbolic network. 

For example, one shore station acting as a master for a hyperbolic 
line of position may also generate a circular line of position. A sec- 
ond station acts only as the slave for the hyperbolic pair. 

Another type of composite system combines elliptical and hyper- 
bolic lines of position as shown in Figure 2.12. A principal advantage 
of this combination is that right angle intersections between families 
of lines are present throughout the network. 

2,4,4 Position Accuracy . 

The position accuracy obtained with any electronic system may be 
thought of as composed of two separate parts: repeatability and pre- 
dictability. Repeatability is defined as the reliability with which 
the electronic method allows the survey vessel to return to a particu- 
lar location on the surface of the ocean specified only by the lines of 
position generated by the electronic system. Predictability is the 
measure of the reliability with which the electronic method defines the 
location of a given position in terms of geographic latitude and longi- 
tude rather than the electronic hyperbolic or range values. 

Repeatability error may be considered as being composed of those 
errors which are present when the system is operated in an idealized 
situation. Components of repeatability are dependent only upon the 
characteristics of the equipment hardware and the precision with which 
it can be operated. A principal component of repeatability error is 
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the uncertainty associated with geometric properties of the system net- 
work, such as line spacing and curve intersection angles. Also includ- 
ed are uncertainties caused by instrument construction characteristics 
which limit equipment measurement accuracy. Another portion of repeat- 
ability error arises from the inaccuracy with which an operator inter- 
prets instrument results. 

Predictability error results when environmental variations and 
changes which occur with the passage of time are imposed on the system. 
Such uncertainties arise when the system is operated in a real situa- 
tion. The principal causes of predictability error are variations in 
propagation velocity and atmospheric refraction. 

Both repeatability and predictability errors are composed of sys- 
tematic and random influences. Repeatability contains a substantial 
amount of systematic variation. Therefore, with careful equipment cal- 
ibration and analysis of the system, repeatability errors can be fairly 
well determined and allowed for. Such is not, however, the case at the 
present time with predictability errors. They are difficult to isolate 
and remove. Empirical formulas have been devised for use in attempting 
to correct for atmospheric condition variations. However, the useful- 
ness of such formulas remains questionable(lO) . 

As mentioned above, repeatability error can be satisfactorily de- 
termined because it depends primarily upon known factors. The develop- 
ment of the expression for the repeatability error of a position which 
follows, applies equally well to a ranging, hyperbolic or composite 
positioning system. In Figure 2.13, S is the apparent location of an 
electronically determined ship position in the geometric network of the 
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system. The angle of intersection of the lines of position, be they 
hyperbolas or circles or one of each, is designated as y, The standard 
errors of the displacement of the two lines of position, denoted as RED 
and GREEN, are m D and m„. 

True curves of constant probability would be ellipses centered at 
S as shown in Figure 2,14. Since for most systems these curves are 
somewhat difficult to compute and construct, for this general develop- 
ment a more suitable figure, the circle, is substituted. The actual 
position of the fix may then be said to lie anywhere within the circle 
centered at S with radius mg with a 68% probability. In this simpli- 
fied form, the radius m g of the error circle is considered to be the 
semi-major ax ■* of the error ellipse. 

The position error of a single observation m^ in terms of m^, m^ 



32 



and the angle y is expressed as follows: 



m: 



m; 



sin 2 y 



m: 
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sin 2 y 



sin 2 y 



(2.9) 



Hence, with a total of n observations: 
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n 



( 2 . 10 ) 



Now if q^ and are the standard errors of each line of position 
obtained by an error analysis of the system concerned, then: 
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Next, a correlation factor R which is one if every radius vector is 
common to at least two pairs and zero if none is used in more than one 
pair is introduced as: 
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Substituting (2.9) into (2.10) gives: 

" 2 - i Zm2 ♦ Ejt| 2 ♦ 2Zm 

(in 2 v L 



mi 



R 



G m R cos Y J 



(2.13) 
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Solving (2.11) and (2.12) for the displacements m^ and m^ and substi- 
tuting into (2.13) yields: 
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(2.14) 



Simplifying (2,14), the standard error of a position, the radius of the 
error circle in Figure 2.13, is: 



m,. 
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(2.15) 
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Finally, the weight to be assigned to each position would be computed 



as: 



(2.16) 



In the case of a hyperbolic system, the error components in the 
directions of the rectangular coordinate axes can easily be obtained. 
Instead of substituting the standard error circle for the ellipse, the 
erroneous displacements of the point S along the RED and GREEN position 
lines are computed. These errors, denoted as e^ and e^ respectively in 




Figure 2,14. Repeatability error ellipse for electronic position- 
ing systems. 



Figure 2.14, are the conjugate axes of the error ellipse. They are 
computed from the geometry of the figure as: 



m 



R 



and 



R sin y 

It can also be shown, according to LAURILA(8), that: 
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(2.17) 
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1 



and 



1 



(2.18) 
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m G = q G J G 



sin S G 

T” 



where 1^ and are the lanewidths on the RED and GREEN baselines and 



q R and q^, are the electronic errors in the RED and GREEN position 
curves as described earlier in the development. Both 1' and q values 
are constant for a given system. 

Now if the error components e R and e^, are considered to be mu- 
tually independent, the following expressions may be written for their 
projections in the X and Y directions where a R and are the angles 
between the RED and GREEN lines respectively and the Y-axis of the car- 
tesian grid system on which the electronic lattice is superimposed. 



■x "V (e R 5i " “g )2 * (e G sln °R )2 

V (2.19) 

a *\J ( e R cos a G ) 2 + ( e G cos « R ) 2 

Combining (2.17), (2.18) and (2.19) and simplifying the resulting ex- 
pressions, the standard errors in the X and Y directions become: 

"x ■ cosec x/(<i r 1r sin v 2 * ( Vg sin v 2 

V (2.20) 

my = cosec ^(q R l R cos a G ) 2 ♦ (q G l G cos a R ) 2 

This development assumes that the error components for the two 
lines of position are mutually independent. No information could be 
found either from a review of the literature or in discussions with 
personnel at the U.S, Naval Oceanographic Office concerning estimation 
of the dependence between survey accuracy determinations at intersec- 
tions of lines of two or more families of circular or hyperbolic curves. 
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Allowing for the future possibility that estimates of covariance 
may be determined, the program for the numerical example calculation in 
Section 4.6 has been designed to accept and utilize such information. 
The following table, extracted from BIGELOW(IO), is intended to 



Table 1, Electronic Positioning System Accuracies 



System 


Type 


Range 

(naut.mi.) 


Accuracy (1) 
(feet) 


Ac curacy (2) 
(feet) 


Loran-C 


hyperbolic 


1200 


50 


1200 


Lorac-A 


hyperbolic 


200 


15 


400 


Decca Survey 


hyperbolic 


200 


25 


300 


Hi-Fix 


hyperbolic 


40 


25 


150 


Lambda 


circular 


425 


25 


250 


DM Raydist 


circular 


200 


12 


100 


Hi-Fix 


circular 


30 


12 


100 


Raydist 


composite 


75 


15 


150 


(1) For best part of service area under excellent conditions. 

(2) At column three range under normal conditions. 



give an indication of the accuracies (repeatability) which can be ex- 
pected with the electronic positioning systems now in use for hydro- 
graphic surveying. These systems are the ones deemed most suitable for 
employment with the application developed in Chapter 4. Of signifi- 
cance is the generally greater loss of accuracy in the hyperbolic 
systems compared with the direct ranging systems at the maximum useful 
ranges. 



2.4.5 Offshore Transmitter Platforms . 

In conventional surveys, station sites for electronic positioning 
system transmitters are placed at known points ashore. These shore 
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points are located at geodetic network stations or are connected to 
them by triangulation, trilateration or electronic traversing. 

In the problem dealt with in Chapter 4, however, the use of land- 
based station sites is precluded in the vicinity of the offshore sur- 
vey. This means that two alternatives are available in so far as the 
use of an electronic positioning system is concerned. First, a long 
range system such as Loran-C could be utilized with a corresponding 
decrease in accuracy as compared with the short range systems normally 
employed in coastal hydrographic survey work. The second choice would 
consist of establishing the transmitters on floating platforms anchored 
offshore, whose positions in turn could be accurately determined. The 
latter of these two alternatives is treated in this discussion. 

The U.S. Naval Oceanographic Office has developed such a technique 
and has demonstrated its feasibility for the operation of short range 
electronic systems along inaccessible coastal regions(ll). For both 
the direct ranging systems with two stations and hyperbolic methods 
with three, the use of floating platforms has been demonstrated suc- 
cessfully. In the former case the platforms are located to provide 
baseline distances of up to twenty-five miles and in the latter lines 
of up to sixty or more miles. In each case, of course, transmitter 
separation distances are dependent upon the specific positioning system 
being employed and the environment of the survey area. 

For such an offshore installation the system transmitters are lo- 
cated in small boats, barges or buoys positioned five to ten miles from 
shore. When boats and barges are used, they are anchored bow and stem 
to provide for maximum positioning stability. Also, to ensure that the 
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floating platform is maintained in its proper position during the 
period of the survey, reference buoys designed to permit only minimum 
horizontal movement are located and emplanted near the floating plat- 
form for position recovery during the survey operation. 

The offshore platforms are positioned by ranging or trilateration 
methods from known positions. Of particular interest to this investi- 
gation is the method used when direct occupation of coastal points is 
impossible. In this situation a helicopter with an airborne position- 
ing system receiver hovers directly over or as near the shore point as 
practicable. Range readings are taken and are repeated by making a 
number of passes over the point to allow computation of as accurate a 
distance as possible when the individual readings are later averaged. 
Another means employed, when hovering is impossible, is to fly over the 
point in a fixed wing aircraft at various widely separated headings, 
preferably about 90° apart. With a minimum of four such passes, on 
each of which the range value is recorded when passing over the known 
point, acceptable mean range values can be obtained. 



2.5 OTHER POSITIONING SYSTEMS . 

In addition to the visual and electronic methods, there are a num- 
ber of other systems which can be used for the control of hydrographic 
surveys. Most are still in the experimental stage of development as 
concerns coastal hydrography applications. However, with the rapid 
technological improvements being made in these and related systems, it 
is not inconceivable that in the near future any or all of them may 
prove capable of supplying the necessary position accuracy required of 
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the ship-based problem treated in Chapter 4. For this reason certain 
of these other systems are briefly summarized. The capabilities and 
limitations as concerns the use of these systems in coastal hydrography 
are also presented. These characteristics are included to indicate why 
and how such systems could either improve or possibly degrade the re- 
sults, compared with information obtained by the conventional visual 
and electronic methods now in use. Of the three methods described, the 
one which appears to show the greatest promise as a locating device for 
coastal surveying is the laser system. 

2.5.1 Photogrammetric Methods . 

Panoramic terrestrial photography taken from the survey ship has 
been utilized on a limited basis in an effort to supplement the usual 
coastal hydrographic survey information. Recently an attempt was made 
to integrate such photography of the coast with an automated survey 
system. In this test, camera circuitry caused interference in digi- 
tizing depth data and skips in data recording. Problems were also en- 
countered in processing and analyzing the photography( 12) . This last 
mentioned deficiency in interpretation of the photography has been the 
greatest cause of problems with such systems. Poor contrast in the 
coastal areas with the resultant difficulty in accurately locating the 
shore targets has thus far limited the usefulness of photography in 
yielding the necessary object position information with which to recon- 
struct the survey situation. 

A second photogrammetric method used experimentally to position a 
survey ship at sea involves an application of photogrammetric satellite 
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triangulation. The survey vessel is equipped with a precision satel- 
lite camera mounted so as to isolate it from the ship's motion. A 
gyrocompass is used to align the leveled camera with geographic north. 
The camera thus oriented is then moved to point in the predicted direc- 
tion of the satellite as it passes overhead. 

At least two other such cameras are mounted at known locations 
ashore and are also pointed in the predicted direction of the satellite 
at the time of its passage. The shore cameras are precalibrated with 
synchronized exposures to obtain their orientation. The photographs 
taken of the satellite from the shipboard camera are used to determine 
its orientation directly from the stars in the background(13) . A syn- 
chronized timing system permits all cameras to photograph the satellite 
simultaneously for the purpose of determining the ship position. 

The system is considered capable of achieving positions within 50 
to 100 feet under normal conditions. An obvious disadvantage is that 
good visibility is required. Also, the data processing to obtain the 
position cannot yet be performed on a real-time basis since a portion 
of the data is at the shore stations. The lack of a continuous posi- 
tioning capability, because of the elapsed time between satellite pas- 
sages, also limits its applicability to survey operations. This defi- 
ciency and the possibilities for overcoming it are discussed in Article 
2,5.3. 



2,5.2 Lasers . 

Considerable effort has been expended in the development of a 
laser theodolite for precise ship positioning in coastal hydrographic 



40 



survey operations ( 14 ) . At the present time such laser systems are 
still in the experimental stage. However, they show considerable 
promise of being adopted on an operational basis with additional re- 
search and improvement. 

The laser systems are to be used for both of the usual visual 
positioning methods, resection from the ship and intersection from 
shore stations. The relative merits of the two alternative methods 
are essentially those previously described in conjunction with the 
conventional visual systems. The primary advantage of the ship-based 
system is the single location of the entire positioning operation. 

The shore-based problem relies for its greatest strength on the stable 
platform and the precise location of the known terrestrial stations 
from which the observations are made. For these reasons the intersec- 
tion solution is again more accurate than the resection from the ship. 

It appears that the ship system could be readily adapted to the 
application described in this investigation since, in a hostile envi- 
ronment, a landmark with good reflectivity may be used as the shore 
target in place of the usual retrodirective target. For the shore- 
based problem, observations are taken from known stations to targets 
mounted on the ship. 

With both systems, range as well as the usual direction informa- 
tion is obtained. With relative bearings to a tenth of a degree and 
line of sight ranges accurate to ten feet, the laser method could con- 
siderably strengthen solutions obtained through the adjustment compu- 
tations developed in the Chapter 4 problem. 

Although the laser beam is affected by atmospheric refraction, the 
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use of dispersion techniques and repetition of readings can substan- 
tially reduce the uncertainty of observational measurements. Compared 
with conventional visual systems, the laser method has the added advan- 
tage of being able to operate in fog, smoke and haze; conditions which 
render the usual visual methods impossible. 

2.5.3 Satellite Systems . 

Doppler satellite navigation systems are now capable of deter- 
mining reliable positions at sea with an accuracy of less than 300 
yards. Utilizing two receivers, the accuracy with these systems can 
be improved to less than 30 yards(l5). To obtain a position, the sat- 
ellite must be observed over a twelve minute period during its passage. 

The lack of a continuous position determination capability is a 
distinct disadvantage of satellite systems as applied to hydrographic 
operations. Intervals of almost two hours, the period of an orbit, to 
more than three hours when the satellite passes within about ten de- 
grees of the observers zenith, are inherent with such a system(16). 

In order to obtain positions on a continuous basis another accu- 
rate system, such as an inertial navigation method capable of giving 
continuous velocity and heading of the ship, must be integrated with 
the satellite method. While satellite systems will provide satisfac- 
tory results for certain types of ocean survey operations, their appli- 
cation to coastal hydrography is considered minimal without significant 
improvements in overcoming the limitations enumerated above. 
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3. AUTOMATION IN HYDROGRAPHY 



3.1 GENERAL . 

The large-scale use of computers and automated procedures in 
hydrographic surveying operations is now being implemented by organi- 
zations such as the U.S. Naval Oceanographic Office and the U.S. Coast 
and Geodetic Survey. Hydrographic survey ships are being equipped with 
the necessary systems and equipment to allow depth, gravity, magnetic 
and positioning data to be acquired and processed automatically. 

With the use of on-line computers, it will be possible to process 
the survey data and produce smooth sheets containing the hydrographic 
data on a real-time basis shortly after the data has been collected. 
From the smooth sheets, field charts will be produced within a matter 
of hours. This feature will be especially useful in areas of military 
operations where charts are required in the minimum length of time. 

Since the use of computers aboard survey ships for the collection 
and processing of hydrographic data is now becoming a reality, the 
adjustment procedures developed in this investigation have been de- 
signed for such processing techniques and are programmed in the For- 
tran IV programming language. 

3.2 AN AUTOMATED SURVEY SYSTEM. 

In order to more clearly illustrate the type of an automated sur- 
vey system to which the adjustment procedures developed in this inves- 
tigation could be applied, the HYDRA System to be employed by the U.S. 
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Naval Oceanographic Office is described briefly. For a complete ex- 
planation of the system, see SPINNING et al (12). In Section 3.3, mod- 
ifications to this system which would be required in order to permit 
acquisition and processing of the data required to solve the Chapter 4 
problem are discussed. 

The HYDRA System was developed principally for military area sur- 
vey applications, also the primary area intended for the ship-centered 
adjustment procedure developed in Chapter 4. As such, this particular 
system is considered to be uniquely suitable as the method by which the 
problem could be implemented. 

3,2.1 Data Acquisition . 

The HYDRA Survey System basically performs three distinct func- 
tions: position fixing, depth acquisition and data recording. For 
position fixing, an electronic positioning system is employed. A digi- 
tizer unit is utilized to translate the output from the shipboard re- 
ceiver into computer language, i.e, to convert lane values into digi- 
tized form. 

A digital depth finding system consisting of an echo sounder, 
transducers and a depth digitizer is used to acquire sounding data. 

As with the positioning system, the digitizer forms the link between 
depth acquisition and data recording. The digitizer converts the tra- 
veling time between transmission pulse and bottom echo, measured by 
the echo sounder, into a digital output value. 

The digitized position and depth data is correlated and sequenced 
by a digital control unit before being recorded on magnetic tape. The 
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information in both digitizers is interrogated and stored together with 
an internally generated time signal. The information is then sequen- 
tially scanned and recorded on magnetic tape in its proper location. 

3.2.2 Data Processing . 

The HYDRA Automated Data Processing System performs two separate 
functions: data processing and plotting. First, analyses are performed 
on the data collected. Errors in the data and weak areas such as loca- 
tions where insufficient information has been collected are determined 
in the edit procedure. The second step entails plotting the ship's 
track by means of the digital plotter in combination with the data 
acquisition system. The third part of the process consists of correla- 
ting the depth information with the positioning data and analyzing the 
results to produce an optimum smooth sheet plot of the sounding data, 
to include the plotting of depth contours if desired. A paper tape is 
generated by the processor for use in plotting the results. This 
plotter tape, when fed into an off-line digital plotter, causes a 
smooth sheet containing the hydrographic data to be produced. 

3.3 APPLICATIONS AND MODIFICATIONS . 

The HYDRA Survey System was developed to satisfy the requirements 
for river surveys in military warfare regions. With modifications to 
allow the acceptance of positioning information from various other lo- 
cation determining systems such as satellite and visual, the use of 
this automated procedure could be extended to other areas of coastal 
hydrographic surveying and also to ocean surveying applications. 
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Two basic modifications would be necessary for the HYDRA System 
to render it capable of incorporating the procedures of the first prob- 
lem dealt with in this investigation, the self-contained system aboard 
the survey ship. First, provision would have to be made for reading in 
the visual data; sextant angles and gyro compass bearings; and corre- 
lating it with the depth and electronic positioning information. Pro- 
cedures have already been proposed for the integration of gyro informa- 
tion into automated acquisition procedures ( 17) . In order to incor- 
porate sextant angle information, a punched card procedure could easily 
be developed to introduce the angular and associated time information 
to the digital control unit for integration and sequencing with the 
other input data. 

As a second modification, or more properly an addition, a program 
to perform the required adjustment calculations would need to be in- 
corporated in the processor library to enable the required computations 
to be accomplished by the computer. Certain other minor modifications 
would also be necessary in order to convert the input data into a form 
suitable to allow the adjustment calculations to be performed. 
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4. DETERMINATION OF SHORE AND SHIP POSITIONS FROM VISUAL 



AND ELECTRONIC MEASUREMENTS AT SEA 



4.1 GENERAL . 

The approach to be taken in this problem utilizes the principles 
and procedures previously described for the ordinary coastal hydro- 
graphic survey. However, it varies considerably in that portions of 
ordinarily known information and available data are now considered as 
unobtainable for various reasons. With these limitations imposed, the 
remaining tools at the disposal of the hydrographer, including high- 
speed computer processing, are utilized to fill in the missing infor- 
mation and arrive at the desired positional results for the survey. 

With this lack of certain usually available data, the accuracy of 
the ultimately determined positions will necessarily be somewhat less 
than had all of the supporting information been present. However, in 
operational situations it often becomes imperative to sacrifice a cer- 
tain degree of accuracy in order to allow the required mission to be 
accomplished in an expeditious manner. 

The situation to be dealt with in this instance is best applied 
along a stretch of inaccessible coastline. In this area a requirement 
exists for a rapid and fairly accurate determination of the positions 
of prominent objects ashore and water depths along the coast. As an 
example of such a situation the following military operation and the 
associated hydrographic requirements are outlined. Water depth infor- 
mation and landmark locations are needed along the coast of a hostile 



47 



nation so that combatant naval vessels may navigate safely and fix 
their positions visually in order to provide accurate offshore gunfire 
support for destruction of enemy fortifications. Accurate inshore 
soundings are also required to facilitate and insure the navigational 
safety of anticipated amphibious landings. 

The only available nautical charts of the area are old and incom- 
plete. They do not contain sufficient information to enable the navi- 
gator to fix his position visually and are practically devoid of 
sounding information. They include only a sparse few topographic fea- 
tures which could be used as navigation landmarks. The location of all 
charted information is questionable owing to the poor accuracies of the 
original surveys. An additional restriction imposed in this situation 
is that the area cannot be overflown for low level reconnaissance pho- 
tography for use in planning and executing the survey. 

Adjustment computation procedures will be applied to the observed 
data in order to arrive at the most likely values for the unknown quan- 
tities in the problem. Estimations of accuracies for all adjusted 
values will also be determined. 

With sufficient shore control, the mechanical solution described 
in Article 2.3.2 normally provides adequate positions for the placement 
of sounding data on nautical charts. However, with little or no known 
terrestrial position information this conventional procedure cannot be 
applied. Also, even with satisfactory visual control ashore, the con- 
ventional sextant resection method may not yield results accurate 
enough to fulfill the most stringent offshore positioning requirements. 
For these reasons the adjustment procedure which follows has been 
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developed for application to coastal survey situations. Although the 
use of adjustment computations has seldom been applied in coastal 
hydrography, procedures have been devised for certain aspects of ocean 
surveying(18) and geodetic positioning at sea away from the coast(19). 

A planimetric approach is taken to this problem since all visual 
observations, both sextant angles and compass bearings, are normally 
measured between objects in approximately the same horizontal plane 
as the observer aboard the survey vessel. When it becomes necessary 
to measure a sextant angle between two points with a significant dif- 
ference in elevation, the angle is corrected before being used in plot- 
ting or computations by a graphical method based on the formula: 

cos 0 



cos a = 



(4.1) 



cos h 

In this expression a is the horizontal or computed angle, 0 is the ob- 
served inclined angle and h is the angular elevation of the elevated 
object from the point of observation. The angular elevation h is nor- 
mally determined by a sextant observation from the survey vessel. 



4.2 PROBLEM OBSERVATION DATA . 

It is assumed that an electronic positioning system has been es- 
tablished so as to give suitable coverage in the coastal area to be 
surveyed. Transmitter sites for the electronic system are established 
outside the inaccessible area, either in neighboring friendly territory 
or, if this is not feasible, aboard floating platforms moored at off- 
shore locations as described in Article 2,4.5. 

The survey ship, as it proceeds along the coast at a distance of 
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two to five miles from shore, obtains the usual sextant angles to prom- 
inent terrestrial objects. If the shore points are at elevations con- 
siderably above the height of the observer aboard the ship, the angles 
are reduced with (4.1) before being entered into the adjustment compu- 
tation, This preliminary reduction of inclined angles, when necessary, 
allows the problem to be dealt with planimetrically in all cases. 

Each time sextant observations are made, a gyro compass bearing is 
taken to the extreme object in the direction closest to the ship's 
track. This particular point is chosen since, in most instances, its 
bearing from the ship will be changing at a slower rate of speed than 
will bearings of the other shore points. A principal reason that only 
one such direction is observed for each ship fix position is that nor- 
mally only a single gyro compass repeater is available for making such 
a reading. 

As the visual data is being obtained, the electronic positioning 
system coordinates of the ship are recorded at each designated fix 
time. Depth information is recorded continuously during the operation. 
Any other desired information such as gravimetric determinations are 
also made simultaneously with the visual observations. Times are noted 
for all data readings to permit later correlation of the information. 

Referring to Figure 4.1, the points through are successive 
positions of the survey vessel proceeding along the coast, The points 
T, U and V represent prominent objects along the shoreline such as 
towers, buildings and natural features between which the sextant angles 
and are observed. All compass bearings t^ are taken to 



50 




4,3 ADJUSTNENT PROCEDURE . 

The generalized least squares adjustment method with parameters 
treated as observables is utilized for the solution of the problem. 
The matrix notations and adjustment procedure are according to 
U0TILA(20) . Owing to the linear relationships which exist among the 
observed quantities and the unknown parameters, and the method used 
for computing approximations for the parameters, a satisfactory solu- 
tion should normally be achieved without the necessity of recycling 
the adjustment. For completeness of the derivation, however, at the 
end of the general adjustment procedure described in this section and 
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in the detailed equation formations of Section 4.4, the modifications 
necessary for recycling the adjustment are explained. 

Sextant angles and compass bearings taken on the ship compose one 
set of observations, L^. These observations are considered as belong- 
ing to the first mathematical structure F^. The adjusted values of the 
observable quantities are designated Lj and the adjusted values of the 
unknown parameters, the ship and shore positions, are X a . There is 
only one observable quantity, a sextant angle or a compass direction, 
in any Fj equation. 

A second set of observations, L^, consists of the successive 
positions of the survey ship, taken as belonging to the mathematical 
structure F 2 . In this structure, L 2a are the adjusted values of the 
ship positions, which are also some of the unknown parameters X . 

a 

The two mathematical structures may be written as: 



Pl(L la , X a ) = 0 



F 2 (I-2 a , X a ) - 0 



(4.2) 



Through the Taylor series development, (4.2) can be expressed as fol- 
lows : 



9F. 3F 

itr v i *'sr- x * F i CL v V * 0 

la a 



3F 3 p 

v 2 *sr- x * F 2< L 2 b . V = 0 

a 



(4.3) 



(4.4) 



where and V ^ are residuals vectors, X Q denotes the approximations to 
the parameters and X represents the alterations which will be applied 
to the approximations to obtain the adjusted values X a « Using the no- 
tations : 
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B 1 = 317 B 2 * Stf 

1 a 

ap. sf 

A i ’ a!C A 2 1 15T 

a a 

"l ' F l< L lb- V W 2 * F 2< L 2b' V 

the differential forms (4.3) and (4.4) become: 



B 1 V 1 + A x x + w ! = 0 


(4.5) 


b 2 v 2 + a 2 x + w 2 = 0 


(4.6) 



Since there is only one observable quantity in any equation of the form 
Fj, Bj is an identity matrix. Hence (4.5) is reduced to: 

V + A X X + Wj = 0 (4.7) 

Further, VI ^ for the first cycle of the adjustment will be a null vector 
because 1^, the observed values of the ship coordinates, and the 
approximate values of the ship coordinates, are the same. Also, since 
there is only one observed quantity, an x or y coordinate of a ship 
position in any F 2 equation, B2 is an identity matrix. Therefore (4.6) 
reduces to: 

V 2 + A 2 X = 0 (4.8) 

The function which must be minimized in order to fulfill the prin- 
ciple of least squares requirement is thus, using the Lagrange multi- 
pliers method: 

$ = V 1 P 1 V 1 + V 2 P 2 V 2 " 2K i (V l + A 1 X + V " 2K 2 (V 2 + A 2 X) ( 4 - 9 ) 
Kj and K 2 are column vectors composed of Lagrange multipliers. is 
the weight matrix for sextant angle and compass direction observations, 
and ? 2 is the weight matrix for the parameters treated as observations. 
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Taking the partial derivatives of <I> with respect to and X and 

setting them equal to zero, the three sets of condition equations which 
must be fulfilled for all observations and parameters are: 



15 Iv7 - Vi - K i ■ 0 
* % * P 2 V 2 - K 2 * 0 



h - -aTk - a"k 2 = 0 

O A 11 2 2 



(4.10) 

(4.11) 

(4.12) 



The expressions (4.7), (4.8) and (4.10) through (4.12) are a sys- 
tem of five sets of equations having five unknown vectors: V^, V^, K^, 
and X, Since the solution of this problem consists of determining 
the X vector, the five set system is first reduced to a system of two 
sets of equations having only two unknown vectors. 

The procedure for this reduction is begun by first solving (4.11) 
for V 2 and substituting the resulting expression into (4.8) to give: 

P~ 2 h 2 * A 2 X = 0 (4.13) 

Next, (4,10) is solved for Vj and the result substituted into (4.7) 
eiVing: P" 1 ^ ♦ AjX * -Wj (4.14) 



. T 

bach term in (4.14) is then multiplied by A^P^ to obtain: 

A 1 K 1 + A i P i A i X = ” A i P i W i (4.15) 

T 

Finally, (4,12) is solved for A^ and the result inserted into (4.15) 
yielding: 

A 1 P 1 A 1 X " A 2 K 2 “ “ A 1 P 1 W 1 (4.16) 



Thus (4,13) and (4.16) are the desired system of two sets of equations 
with unknown vectors K 2 and X. 
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T 

In (4.16) the term is the normal equation coefficients 

T 

array which will be designated Nj. Also, the term AjPjWj is the con- 
stant vector of the normal equations, to be designated U. Using these 
notations, (4.13) and (4.16) become: 

P" 1 ^ + a 2 X = 0 (4.17) 

NjX - = -U (4.18) 

Solving (4.17) for K 2 , 

K 2 = -P 2 A 2 X (4.19) 

which, when substituted into (4.18) gives: 

N 1 X * A 2 P 2 A 2 X = " U (4.20) 

Factoring out X, 

(N 1 + A 2 P 2 A 2 )X = ' U (4.21) 

But each F 2 equation contains only one unknown parameter, an x or 

y ship position coordinate. Hence the only non-zero elements of A 2 are 

- l’s. Further, by virtue of the F 2 structure, A 2 is a square array 

with non-zero values occurring only in the main diagonal elements 

corresponding to ship parameters (see Article 4.4.12). Therefore 
T 

A 2 P 2 A 2 * s merel y P 2 non-zero elements only on the main diagonal at 

positions where dependence is estimated between ship position coordi- 
nates as explained in Article 4.4.14. In matrix notation then: 

A^P 2 A 2 . P 2 C4.22) 

Substituting (4,22) into (4.21), the factor (N^ + P 2 ) in the resultant 
expression, known as the generalized normals matrix, is denoted as N^ en . 
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Using this notation, (4.21) can be rewritten as: 

(N, ♦ P.JX * N X = -U (4.24) 

v 1 2 J gen v 

The approximate values X Q of the shore point coordinates and the 
ship station coordinates are then combined with the alterations deter- 
mined from (4.24) to obtain the adjusted values of both the terrestrial 
points and the ship station coordinates. 

X = X + X (4.25) 

a o v ' 

If the magnitudes of the final alterations determined by (4.24) 
indicate that the first approximations X q are not sufficiently accu- 
rate, the adjustment would be recycled. In this situation, the X 

3L 

computed from (4.25) for the first cycle would become the X q for the 
second cycle of the adjustment. The adjustment would then proceed as 
outlined above with certain modifications. 

The misclosures vector would no longer be null but would in- 
stead consist of the differences between the initially computed adjust- 
ed x and y ship coordinates and the originally observed values. With 
this addition, the W 2 vector would be included in equations (4.8), 

(4.9) and (4.13). Thus the constant vector U would become: 

U = AjPjWj - P 2 W 2 (4.26) 

With the value of U from (4.26) substituted into (4.23) the problem is 
solved as previously explained to obtain the new alterations X. In 
(4,25) the finally adjusted parameters computed from the second cycle 
would consist of the parameters X plus the approximations X q which are 
now the adjusted values computed from the first cycle. 

Should an additional recycling be deemed necessary, the adjusted 
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values from the second cycle would be utilized as the approximations 
for the third cycle and the operation repeated as described above. 

4.4 FORMATION OF EQUATIONS . 

The quantities necessary to accomplish the adjustment presented in 
the previous section are now described, and the expressions which are 
required in order to perform the indicated computations are developed. 

The network parameter notations to be employed are as shown in 
Figure 4.2. The cartesian coordinates of each ship position are 
denoted as (x^, y^) . For each shore point the coordinates correspond 
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to the point concerned. For example , for terrestrial object V the 
coordinates are (x v , y^) . For the purposes of this general development 
is is assumed that there are three terrestrial points and n ship posi- 
tions. 



4.4.1 Observed Quantities - First Structure . 

The quantities are the directly observed sextant angles 

ai b* B lb» “ 2 b» B2 b* * * * “"b* Bn b 
and the directly observed compass bearings 

t lb» t 2b* ‘ ‘ ‘ l nb 

to the terrestrial point closest to the direction of ships travel. 

4.4.2 Observed Quantities - Second Structure . 

The quantities are the ship position coordinates as determined 
by an electronic positioning system. 

( x lb» n b )»( X 2 b» » • • • C^b* 

4.4.3 Approximate Values of Parameters - First Structure . 

These parameters are the approximate coordinates of the shore 

points. They are computed in the following manner. With the observed 
coordinates of ship positions and S^, the coordinates of point T 
are determined to illustrate the procedure (See Figure 4.2). Using 
plane geometric and trigonometric relationships, the approximations 
for tne distances between ship positions and for the angles P and y 
are calculated. With these values, the approximate distances from the 
snip positions to T are computed according to the law of sines. The 
distance from Sj to T is thus given by: 
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si " P 1T K D ln 



(4.27) 



Dist 



IT 



sin y 



In 



With the observed compass bearing t^y to terrestrial point V and 
the observed sextant angles and 8j, the approximate azimuth t^ is 
computed , 



t lT * *1V “ (a l + 



(4.28) 



The horizontal and vertical increments Ax^ and Ay^ are then written 



as: 



Ax 1T = Dist 1T x sin t 1T 



Ay 1T = Dist 1T x cos * 1T 

The approximate coordinates of point T are thus expressed as: 

*T 0 - x l 0 * Ax lT 
y T 0 * * ^lT 

In similar manner the approximate coordinates of points U and V are 
determined. 



(4.29) 



(4.30) 



4.4.4 Approximate Values of Parameters - Second Structure . 

For ship coordinates, these parameters are the observed quantities 
described in Article 4.4.2. 



4.4.5 Computed Values of sextant Angles . 

The general expression for a computed sextant angle formed between 



shore ncints J and K, subtended at ship position is given by: 



= tan 



X J - x i 


„ -1 
. tan 


X K “ X i 


/j - y i_ 




Jk ' y i 



(4.31) 
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4.4.6 Computed Values of Compass Bearings . 

The general expression for a computed compass bearing from ship 



position to shore point J is: 



t. 

1 



tan 



-1 



X J ‘ x i 



y J * y i 



(4.32) 



4.4.7 First Mathematical Structure . 

This is the model of the form F-CLj , X ) = 0. An expression is 

A 3> d> 

developed for each sextant angle and compass bearing. 

-1 y ' u a ~ '' 1 a ^ -1 ” x l a 



tan 



tan 



-1 



tan 



tan 



-1 



x Ua 


“ x la 


y y a 


y la_ 


> 

X 


‘ x la 


y V a 

• 


1 

M 
j 03 




- *n 

"a 


y Ua 


- y na _ 


x v a 





- tan 



- tan 



-1 




" °l s 



- 61. 



ss 0 



• • 



yv a - y n 



a 



- tan 



tan 



tan 



-1 



tan 



-1 



• 

X T a " x n a 


• 

_ 


y T a ‘ y n a 


“n; 


^a • X n a 


- s 


yu a - yn a 


p n 


x V a “ x l a 


- 1 <• 


y v a " y i a 
• 


• 


• 

• “ 
x v a - x n a 


• 

• 

_ t 


y V a * y n a 





• • 



‘a 



= 0 



(4.33) 



4.4.3 Second Mathematical Structure. 

rh-us model, the parameters treated as observations, is expressed 

as F,(L2 , X ) = 0. The shore point parameters are only treated as 
a a 
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observations in the sense of being observations with zero weight. Thus 
shore coordinate functions do not enter this structure. To present the 
ship parameter portion of the structure in a more meaningful manner, 
the relationships L2 a * ♦ V ^ and (4.25) are substituted permitting 

the observation equations to be written. They take the following form: 

x. ♦ v, - Ax. » 0 

*b 1 

y, + v 2 - Ay « 0 

b *2 1 

x 2. * v 2_ - 4x 2 * 0 
D 3 

>'2, * V 2 - * 0 

b 4 



(4.34) 



x ♦ v, - Ax » 0 

"b 2 2„-l n 

V * V 2, - % ■ 0 

b 2n 



4.4.9 Misclosure Vector - First Structure . 

The misclosures are obtained by evaluating the mathematical 
model with and X q . To illustrate the procedure, the expressions 
for the first sextant angle and the first compass direction are shown 
below. Since there are three observations at each of the n ship posi- 
tions, Wj consists of a total of 3n elements. 



o’ tan 



Xm - Xi 

l) o A o 



yu c - yi 0 



- p* tan" 



tan 




“lb 



= w. 



(4.35) 



-*»;-«i2n.l (4 - 36) 
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4.4.10 Misclosure Vector - Second Structure. 



As explained in Section 4.3, the W 2 matrix is null for the first 
cycle since L2^ and X q are the same quantities. However, for the re- 
cycled adjustment for ship parameters, (xj o , y^ Q ) are the adjusted 
coordinates computed from the first cycle and (x^* yi^) are the °* J " 
served values from the electronic positioning system. In this case, 
that portion of W 2 corresponding to ship coordinates takes the follow- 
ing form: 



x 



1 

o 





y 



i 

o 




X 



2 



0 








(4.37) 





w. 



w. 



2n-l 

2n 



4.4.11 Direction Coefficients - First Structure . 

This array, of the general form A, = 3F../3X , is obtained by 

1 13 

taking the partial derivatives of the structure expressions with res- 
pect to each of the unknown parameters. For ease in computation and 

o 

presentation, the a array is partitioned into two parts, A represent- 
ing the partials with respect to the terrestrial object parameters and 
A designating the partials with respect to ship position parameters. 



A 



1 




(4.38) 
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Developing first the A array, the partial derivatives for a^ and 



tj are presented to illustrate the types of expressions which make up 
the elements of the matrix. 

- pM y T - yp 



3a i 

9xT 



(Xj - Xj ) 2 + (y T - y x ) : 

+ p'(x t - Xj) 



9a 






( X T - Xj) 2 ♦ (y T 

+ p'fru - y x ) 



y l> 2 



(x u ‘ X 1 )2 + (y u - y i ): 



9a. 



- P'CXy - Xj) 



r u 



9t 



1 

*V 



(Xu - x,) 2 . (y„ - yj) 2 

- p'(y v - yp 

(Xy - Xj ) 2 ♦ (y v - yp 2 



9t 



1 



+ p'(x. 



x l ) 



(X, 



x l )2 * (y V " y l )2 



( 4 . 39 ) 



Since partials are taken with respect to the x and y coordinates 

o 

of each of three shore points, the A matrix has a column dimension of 

six. Because there are three observations, two sextant angles and a 

o 

compass bearing, at each of the n positions of the ship, the A matrix 
has a row dimension of 3n. The arrangement of the partial derivative 
elements is shown in the schematic array on page 64. 
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oooo lin. 

Next, the A array is developed. Again the partial derivatives for 
o,j and tj are presented to illustrate the expressions of which the 
matrix is composed. 

9a i p'Cyu “ yp p'Cy T - y x ) 

3x l (x U * X l )2 * (y u “ y l )2 (X T ‘ x i )2 * (y T “ y l )2 

SOj P'CXy - Xj) P’(X T - Xj) ( 4 - 41 ) 

3y l (X U ’ X l )2 ■* (y U ' y l )2 (X T ’ X l )2 + ^ y T ‘ y l ^ 2 
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p'(y v - /j) 



at. 



ax. 



at. 



(X x/ - X,) 2 4 (y y - yp 2 



V V 

p'(x. 



* ♦ 



Xj) 



(4.41) 



8 y t (* v 



x l> 2 * (y v - ’’l’ 2 



Because partials are taken with respect to the x and y coordinates 
of each of the n ship positions, the A matrix has 2n columns. Since 
there are three observations at each ship position, the row dimension 
is 3n. The matrix array of the partial derivative elements is: 



3n^2n 



3ot l 8ot l o 

■5x7 WT 

ag i ag i 0 

"5x7 "577 0 



• • • 



• • • 



9a2 3oi2 

3x7 ay7 * * * 



o 9B 2 as 2 

3x 2 3y 2 



• • • 



• • • 



• • • • • 



° • ••& & 



Xn ay n 

98n ag n 
3x n ay n 



3ti 3tj « 

177 ° 



3 t2 8t2 

ax 2 ay 2 



• • « 



• • • 

• • • 

• « • 

• • • 



0 atn at n 

0 * * * W 177 



(4.42) 
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4.4,12 Direction Coefficients - Second Structure. 



This array, of the general form A_ * 3F-/3X , is obtained by 
taking the partial derivatives of the structure expressions with re- 
spect to each of the parameters. Since the shore point coordinate 
parameters are treated as observations with zero weight, the partial 
derivatives for these parameters are simply zero. For ship coordi- 
nates, the partial derivatives are: 

3f i 3(-x.) 

= -1 



3Xj 


3x l 


3f 2 


K-yj) 


*1 “ 


3 ^i 


3f 3 
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3(-x 2 ) 


3x 2 


* X 2 


Sf 4 


3(-y 2 ) 
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* 

• 


• 

• 

3f 2n-l 


• 

• 

3 (-X ) 
v n' 


&X 


ix 


n 


n 


3f 2n 


3(-y„) 




3 ^n 



(4.43) 



-1 



= -1 



Since tho partial derivatives are taken with respect to the x and 
y coordinates of each of three shore points and n ship positions, the 
matrix has a column dimension of 2n*6. Because there are two obser- 
vations, an x and y coordinate, at each of the three shore points and 
n ship positions, the A 2 matrix also has a row dimension of 2n*6. Fur- 
ther, it is a diagonal array with the only non-zero elements (-l's) oc- 
curring in positions corresponding to ship coordinates. 
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4.4.13 Weight Matrix - First Structure . 

For sextant angles, it is assumed that the standard errors for all 



observations are equal . 

directions. That is, 

W * m a 
°1 Si 



The same assumption is also made for compass 



» m * m fl . . 
02 D2 

= m. . . . » m. 
t 2 t 



m 

a 



(4.44) 



Since all observations are considered to be independent, Pj is a diag- 
onal matrix of the reciprocals of the sextant angle and compass bearing 
standard errors squared. Since there are three observations at each 
ship position, the square dimension of P^ is 3n. 



4.4.15 Weight Matrix - Second Structure . 

?2 is the weight matrix for all observed parameters. Since the 
shore coordinates are treated as completely unknown parameters (obser- 
vations with zero weight), ? 2 takes the general matrix form: 



2n*6 2n*6 



X 



-1 



(4.45) 



£- is the variance-covariance matrix associated with the ship coordi- 
nates. In its most generalized form, it would be entirely filled with 
variance and covariance estimates. If, however, it is assumed that 
systematic errors which would affect all x and y ship coordinates can 
be ascertained and removed or that the positioning method is being 
utilized as a relative "localized" system, the successive x coordinates 
and successive y coordinates can be considered as independently deter- 
mined quantities. The array would then consist of "n" 2x2 blocks on 
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the diagonal, each block representing a ship position. The four ele- 
ments of each block are the estimated variances of the x and y coordi- 
nates on the principal diagonal and the covariance between x and y in 
the off diagonal positions. If for each ship position the x and y 
coordinate determinations are considered as independent quantities, 
then E- becomes a diagonal matrix. Further, if all coordinate values 
are assumed to be of equal accuracy, is a unit matrix. For the 
localized system 2x2 block form which is utilized in this problem, E- 
is represented as follows: 
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2n x 2n 
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(4.46) 



4.4.15 Normals Matrix - First Structure . 

.th the Aj matrix partitioned into A and A, the normals matrix N 
is formed as follows: 



T 

N « A P A = 

1 1 11 



OT 

A 


Mr? : il . 


Or O ! OT - 

A* PjA ; ATjA 


a t 


L p ilL A • A J 


-t o ; _t - 
A l PjA : ATjA 



(4.47) 



Now, setting A*PjA equal tc N, the expression for Nj may be written as: 
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(4.48) 



N. 
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4.4.16 Generalized Normals Matrix. 



The generalized normals matrix N g en is formed by adding P 2 to Nj 



The square dimension of N is 2n+6. 

gen 



N «* 
gen 
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zo 
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Zi 
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N . Ej*» 



(4.49) 



4.4.17 Constant Vector of Normal Equations . 

The vector U is formed as explained in Section 4.3 wherein the 
component matrices are as described in previous articles of this sec- 
tion. For both the first and recycled adjustment, U consists of 2n+6 
elements representing the x and y coordinates of each shore point and 
ship position. For the first cycle, U is formed as follows: 



U = AjPjWj * 



OT 

A 



a t 



NM 



°T 

*Vi 



*Vi 



(4. SO) 



When the adjustment is recycled, U becomes: 



U = A^PjWj - P 2 W 2 = 



OT 

A Vi 



A Vi 






(4. SI) 



4.4.18 Adjusted Parameters . 

The corrections X obtained by solving (4.24) are combined with the 
approximate values of the shore points and observed ship positions as 
shown in (4.25’ to obtain the finally adjusted values X of the 

a 
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terrestrial objects and ship position coordinates. Since there are n 
ship positions and three shore points, each with an x and y coordinate. 



the X & vector is composed of 2n*6 elements. 



Xj 

a 




*T 

o 


♦ 


Ax^, 


y T 

a 




y T 

0 


♦ 


Ay T 


x u 

a 




*u 

o 


♦ 


Ax u 


y U 

a 




/ U 

o 


♦ 


Ay u 


X V 




x v 


♦ 


Ax V 


a 




V 

O 




V 


y v 

a 


a 


y v 

O 


♦ 


Ay v 






x i 

0 


♦ 


Ax l 


•'i 

a 

• 




y l 

0 

• 


♦ 


Ay l 

• 


• 

• 

X 

n 

a 




• 

• 

x n 

0 


♦ 


■i 


\ 






♦ 


Ay n 



(4.52) 



4.4.19 Adjusted Sextant Angles and Compass Bearings . 

These values, Lj , are obtained by adding the residuals vector V, 
a x 

computed using (4.7) to the originally observed sextant angles and 

comn'ss directions L^. Since there are three such observations at 

each ship position, consists of 3n elements. 

a 
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1, 




♦ V„ 


1 




1, 1 


a 




b 


l 2 




\ * V 2 


a 




b 


1, 




1, ♦ v_ 


3 




3, 3 


a 




b 




s 




• 




• « 


• 




• a 


• 




• • 


1 3n-l 




1 ♦ V 

3n-l. 3n-l 


a 




b 


1 3n 




V * v 3n 


a 




b 



(4.53) 



4.S ERROR ANALYSIS. 

In order to examine the accuracies of the adjusted values of the 
parameters and the observed quantities, the variance-covariance ma- 
trices for each set of values is obtained. The sequence of operations 
necessary to arrive at these arrays is explained below, continuing the 
matrix notations utilized in the previous two sections. 



4,5,1 Variance of Unit Weight. 

2 

The variance of unit weight, denoted as m Q , is first determined. 
All weighting in the problem is based on the same variance of unit 
weight. This unit weight variance is expressed as the sum of squares 
of the weighted residuals divided by the degrees of freedom. 



T 

V PV 



2 = 

m o ~ degrees of freedom 



(4.54) 



The weighted residuals square sum for the entire problem is the sum of 
the square sums from each of the two mathematical structures. 



V T PV = VjPjVj + V^P^ 



(4.55) 
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In partitioned matrix form, (4,55) is expressed as: 




(4,56) 



Pj is explained in Article 4.4.13 and P 2 in 4,4.14. The square di- 
mensions of Pj and are 3n and 2n+6 respectively. Vj is computed 
from (4.7) and is a column vector of 3n elements. V 2 is calculated 
using (4.8) and consists of 2n+6 elements. The last 2n elements are 
simply the negative X vector of alterations to ship coordinates. 

The number of degrees of freedom is determined as follows: 

degrees of freedom = r ♦ s - u (4.57) 



where 



r * the number of sextant angle and compass observations, 

s ■ the number of ship position x and y coordinates, 

u = the number of shore point x and y coordinates. 

For the general problem with three shore points and n ship positions, 

r » 3n, s » 2n and u * 6 (4.58) 

Hence the number of degrees of freedom, substituting (4.58) into 
(4.57), is 5n - 6. Thus the variance of unit weight is finally ex- 
pressed as: 



T 

VP V ♦ 

2 1 11 

m o 5n - 6 



T 

VP V 
2 2 2 



(4.59) 



4.5.2 Variance-Covariance Matrix of Parameters . 

The weignt coefficient matrix of the unknown parameters is simply 
the inverse of the generalized normals matrix developed in Article 
4.4.16. 
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(4.60) 



N 



gen 



Thus the variance -covariance matrix Z ^ which indicates the accuracy of 



the adjusted values of the parameters is: 



'X = m o Q X 



(4.61) 



It is a square array of dimension 2n+6. 



4.5.3 Variance -Covariance Matrix of Observations . 

The weight coefficient matrix for the adjusted values of the ob- 
served sextant angles and compass directions is: 



% • a i n ; 1a I 

A a 



(4.62) 



where and are as described in Articles 4.4.11 and 4.4.15 respect- 
ively, Its square dimension is 3n since two sextant angles and a com- 
pass direction are observed at each ship position. Multiplying (4.62) 
by the variance of unit weight from (4.59): 

Z. * m 2 AX^T (4.63) 

Li O 1 1 1 

A a 



4.6 NUMERICAL EXAMPLE . 

The unavailability of actual hydrographic survey data in a form 
which would lend itself to adaptation to this problem necessitated the 
compilation of fictitious information with which to illustrate the 
computation procedure. In order to simulate the actual observation 
situation as closely as possible, the following method was employed. 

A 28" x 30" plotting sheet prepared by the Naval Oceanographic Office 
was used as the basis for generating the required information. The 
plotting sheet contained a circular lattice net, longitude and latitude 
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and a UTM X-Y grid coordinate system. Since standard computer program 
routines are available for converting coordinates from any of these 
three systems to the other two, the X-Y grid was utilized for the pur- 
poses of this example. In a practical situation, the electronic 
positions would be recorded and later converted with a computer program 
to the desired X and Y coordinate values. 

The computational procedure for this problem was programmed in the 
Fortran IV language and run on the IBM 7094 computer at the Numerical 
Computation Laboratory of The Ohio State University. Free electronic 
computer services were made available by the University for all pro- 
gramming carried out in this investigation. 

Considerable effort was expended in designing the program for this 
problem in as general a form as possible to allow for variations in the 
types and amount of input information. Many combinations of ship 
headings and observation data were tested to insure that the program 
would produce proper results for all conceivable situations. The nu- 
merical example presented in this section is accomplished with three 
shore points and four ship positions. However, the number of each type 
of position can be increased as desired depending upon the amount of 
input data to be processed. 

4,6.1 Problem Input Information. 

(a) Coordinates for each of four ship positions. 

Source: arbitrarily placed on plotting sheet. 

Data units: meters. 
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Ship Position 


X 


Y 


1 


299865.0 


1960800.0 


2 


299992.0 


1961573.0 


3 


300081.0 


1962378.0 


4 


300224.5 


1963203.0 



(b) Sextant angles and compass directions observed at each ship 
position. 

Source: Measurements made on plotting sheet with three-arm 

protractor. 



Position 


a 


6 


t 


1 


31° 33' 


24° 25* 


32590 


2 


32° 04* 


© 

o 

o 

cr> 

<N 


316?5 


3 


29° 37' 


32° 36* 


305?7 


4 


24° 35* 


32° 20* 


291?2 



(c) Side of survey ship from which angle and direction observa- 
tions are made: Port side. 



(d) Variance and covariance estimates for electronic positioning 
system coordinates of ship positions. 

Source: Arbitrarily selected, based upon information ob- 

tained from various publications describing accuracies achievable with 
short-range electronic positioning systems. 



Data units: meters squared. 



Ship Position 
1 
2 

3 

4 



X variance Y variance Covariance 



15.1 18.6 +2.1 

20.8 19.4 +2.4 

17.9 19.9 +1.8 

16.1 20.8 +2.0 



(e) Standard errors of sextant angle and gyro compass observa- 
tions. 
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Sources: Sextant angle - from estimates given in JEFFERS(7) 

and FAGERHOLM and T1UJNBERG(21) . Gyro compass bearing - arbitrarily 
selected based on discussion presented in B0WDITCH(3) . 

Sextant angle observation standard error: 1 lo 

Gyro compass observation standard error: 5 1 0 



4.6.2 Auxiliary Information Computed by Program . 

(a) The azimuth from each ship position to each succeeding 



position. 



12 


009?330 


*14 


008?509 


*24 


008?118 


13 


007?794 


*23 


006?309 


*34 


009?867 



(1>) Approximate angles y subtended at each shore position between 
each set of ship positions, A partial listing, for the angles between 
ship positions 1 and 4, is as follows: 

At point T At point U At point V 

34?7S0 41?717 33?800 



(c) Approximate coordinates of each 


shore point (in meters) 


Point 


X 


Y 


T 


296810.8 


1960748.4 


U 


296904.6 


1962549.6 


V 


297421.0 


1964290.4 



4.6.3 Adjusted Values and Associated Standard Errors . 



(a) Ship positions (in meters) : 



Position 

1 



X 


m 


Y 


m 


299864.8 


1.6 


1960798.2 


1.6 


299990.8 


1.2 


1961574.4 


1.4 


300082.2 


1.3 


1962380.7 


1.3 


300224.8 


l.S 


1963201.7 


1.6 



2 

3 

4 
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(b) Shore points (in meters): 



Point 


X 


m 


Y 


m 


T 


296810.2 


7.0 


1960747.8 


2.3 


U 


296902.1 


4.1 


1962549.7 


1.9 


V 


297416.6 


4.1 


1964292.4 


3.2 


(c) Sextant 


angles: 








Position 


o 


m 


8 


m 


1 


31° 33 I 0 


i!o 


24° 25 !2 


o!7 


2 


32° 03! 7 


o!9 


29° 00!0 


o!8 


3 


29° 37 '.7 


i!o 


32° 35! 7 


i!o 


4 


24° 34! 6 


i!i 


32° 20!2 


o!9 


(d) Compass 


directions: 








Position 


t 


m 






1 


325?02 


5!o 






2 


316?44 


5!o 






3 


305?75 


5!o 






4 


291?17 


5!o 
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5. DETERMINATION OF SHIP POSITIONS FROM VISUAL 



MEASUREMENTS AT SHORE CONTROL POINTS 



5.1 GENERAL . 

This problem consists of locating the ship position by directions 
observed from the shore, as described in Article 2.3.1. For the pur- 
poses of the adjustment, it is assumed that a series of ship positions 
are to be determined in order to either calibrate an electronic posi- 
tioning system or to determine precise sounding locations for a nauti- 
cal charting or engineering project. The ship will be observed from 
a number of shore stations at specified time intervals as it proceeds 
slowly along the coast in the survey area. 

In practice, directions are usually observed to the ship with a 
continuous tracking azimuth instrument which permits readings to one- 
hundredth of a degree to be made. The instrument is initially oriented 
on another known station. The angle to the ship is then measured to 
obtain the ship direction. This type of instrument is advantageous for 
the situation where the ship is moving at a fairly rapid rate of speed 
as the reading is being taken. It does not, however, permit the maxi- 
mum possible direction accuracy to be achieved. In order that this be 
accomplished, a theodolite capable of considerably greater precision 
would be employed. For the purposes of this problem, where maximum 
accuracy is desired, it is assumed that the theodolite is used with the 
vessel either anchored or moving slowly enough to permit satisfactory 
directions to be observed to it. 
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The ship positions to be determined may be considered as being in- 
dependent on one smother since no connections are made between them. 
Thus each one can be determined by a separate adjustment, utilizing in 
each case all of the information observed at the shore points which 
relates to the ship position concerned. 

5.2 PROBLEM OBSERVATION DATA . 

Referring to Figure 5.1, at each occupied shore station (T, U, V 
and W) two directions are observed to known network stations (1, 2, 3, 
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4 and 5). The directions observed to the known stations , together with 
azimuths computed using the coordinates of the known points, are em- 
ployed to calculate orienting angles for the later determination of the 
orientation of intersecting azimuths to the ship. 

The second part of the problem consists of azimuths being simul- 
taneously observed to the ship target S from the occupied shore sta- 
tions T, U, V and W at time intervals dependent on the ship's speed and 
the survey requirements. Each time a set of directions is observed to 
the ship, the electronic positioning system coordinates of the ship are 
also recorded, as well as depth information and any other hydrographic 
data desired. 

5.3 ADJUSTMENT PROCEDURE. 

The method of variation of parameters is utilized to solve this 
problem. The procedure employed for the intersection adjustment is 
based on a solution given by RICHARDUS(22) . Matrix notations according 
to UOTILA(20) are used to present the adjustment procedure. 

The directions observed from the occupied shore stations to the 
ship compose the observed values of the observable quantities, L^. 

The adjusted values of the observable quantities are designated as L . 
The unknown parameters are the coordinates of the ship position S. The 
approximate values of these coordinates, computed before the adjust- 
ment, are denoted as X . The values obtained as the result of the 

o 

adjustment are X . The amount by which X must be altered to give X 
a o a 

is the correction vector X, 

X » X ♦ X (5.1) 

a o 
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The mathematical structure which is selected for the adjustment 
is a system of equations expressing the adjusted observed quantities 
as functions of the adjusted parameters. 

L a - F(X a ) (5.2) 

Using (5.2) with the computed approximations X q of the unknown para- 
meters, the numerical values of the observed quantities can be rigor- 
ously computed through the mathematical structure. 

L - F(X ) (5.3) 

o v o 

Now the difference between L q from (5.3) and the observed values 
1^ is the corrections vector L. 

L » L q - 1^ (5.4) 

The residuals vector of observations, V, is expressed as the difference 
between the adjusted and observed observations defined above. 

V * L a ’ H* (5,5) 

But it has been shown that the adjusted observations are a func- 
tion of the adjusted unknown parameters (5.2). Therefore, the residu- 
als (5.5) can be written as: 

V - F(X a ) - ^ (5.6) 

Now it was also explained that the adjusted unknowns are the approxi- 
mate parameters plus corrections (5.1). Hence (5.6) can be expressed 

AS • 

Lh ♦ V » F(X o ♦ X) (5.7) 

Linearizing using the Taylor series development: 

^ ♦ V - F(X q ) ♦ AX (5.8) 

where A is the direction coefficients matrix, the partial derivatives 
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of L with respect to X. 

Substituting (5.3) into (5.8) and using (5.4), the result is sim- 
plified to: 

V * AX + L (5.9) 

which is the expression for the observation equations. 

The function which must be minimized in order to fulfill the prin- 
ciple of least squares requirement is now: 

* » V T PV (5.10) 

where P is the weight matrix of the observed quantities. Substituting 
(5,9) into (5,10): 

4> * (AX + L) T P(AX + L) (5.11) 

Expanding the right side: 

$ * X T A T PAX ♦ X T A T PL ♦ L T PAX ♦ L T PL (5.12) 

T 

Now since X and L are vectors and P is symmetric (P = P ) : 

X T A T PL ■ L T PAX (5.13) 

Thus (5.12) can be written as: 

♦ » X T A T PAX ♦ 2X T A T PL ♦ L T PL (5.14) 

T T 

Designating the normals matrix A PA as N and the constant vector A PL 

as U, (5.14) takes the form: 

* * X T NX ♦ 2X T U ♦ L T PL (5.15) 

Taking the partial derivatives of ♦ with respect to X and setting 

T 

them equal to zero, the condition which must be fulfilled for V PV to 
be a minimum is: 

h -jUr * xV ♦ U T • 0 (5.16) 

dX* 

or 

NX ♦ U • 0 (5.17) 
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By solving (5.17) for X, the corrections to the approximate values of 
the unknown parameters are: 

X « -N -1 U (5.18) 

Substituting these corrections into (5,1), the desired adjusted values 
of the unknown parameters, the coordinates of the ship position, are 
obtained. 

The adjusted directions to the ship are obtained by solving (5.9) 

for V and substituting the result into (5.5) to permit the computation 

of L . 
a 



5.4 FORMATION OF EQUATIONS. 

The quantities necessary to accomplish the adjustment procedure 
presented in Section 5,3 are now described and expressions required are 
developed. The notations used are as indicated in Figures 5.2 and 5,3. 



5.4.1 Approximate Values of Unknown Parameters . 

The method employed for the computation of approximate ship coor- 
dinates is similar to the procedure described in Article 4,4.3 except 
that the terrestrial and ship information is interchanged. 

Using plane geometric and trigonometric relationships, approxima- 
tions for the angles labeled Q and y in Figure 5.2 are first deter- 
mined. With the computed distances D between shore positions, the ap- 
proximate distances from shore stations to the ship are computed by 
utilizing the law of sines. For example, from shore point T to the 
ship S: 



Dist, 



sin x 



TS 



(5.19) 



sin y 



TW 
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For the computation of approximate coordinates of the ship posi- 
tion S, the azimuth from terrestrial point T to S is denoted as t T g. 

The increments Ax,pg and are then expressed, using (5.19), as: 

AXjg * Dist Ts * sin t TS 

(5.20) 

Ay TS * Dist^ x cos t^ 

Thus the approximate coordinates of S are: 

*0 ■ X T * “TS 

, (5.21) 

* »TS 

5.4.2 Mathematical Structure . 

The mathematical structure L * F(X ) , the adjusted observations 

3 
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as a function of the adjusted coordinates of the ship position, is now 
formed. It consists of the adjusted computed directions from each 
shore point to the ship position S, 



t TS a tan 
a 



y a " y T 



l US * tan 
a 



x a " X U 



t vs « tan 
a 



*KS ' ta " 
a 




(5.22) 



5.4,3 Approxim ate Values of Observed Quantities , 

By substituting the approximate unknown coordinates from (5.21) 
into the mathematical structure (5,22), the approximate numerical 
values of the observed quantities, L q , are rigorously computed. 



5,4.4 "Observed" Quantities . 

Next the expressions for the observed quantities, L^, are deter- 
mined. It should be noted that as referred to here is not simply 
directly observed directions, but a combination of the observed direc- 
tions and orienting angles computed from the directions to known net- 
work stations. First the orienting angles CK from each occupied shore 
station to two other known positions (see Figure 5.3) are calculated. 
These orienting angles are expressed as shown in (5.23) where the term 
for the computed azimuth from any occupied station i to any known 
network point j is: 



tan 



x. - x. 

-j i 



y i ' y i 
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CL 



(B. 



Tl 



t Tl) * (^T2 



*T2^ 



CB 



U 



U2. 



t U2 ) * (0 U3 



t U3 ) 



(6 



V3 



t V3 ) + (8 V4 



V 



(8 



W4 



*W4^ + C 0 W 5 



t lV5 ) 



W 



(5.23) 



The desired "observed" values are then obtained by adding the 
observed ship directions to the orienting angles computed in (5.23). 
This vector, designated L^, is shown on the following page. 
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h, 



X * ° T 

'us b * °u 

t VS. + °V 

D 

X * ° w 



(5.24) 



5.4,5 Direction Coefficients . 

The expressions for the elements of the direction coefficients 
array A are expressed symbolically as: 



A 



31 

ax 



ii 


a !2 


21 


a 22 


31 


a 32 


41 


a 42 



(5.25) 



where the individual elements, 
following forms: 




illustrated by a^ and a 12 , 

p,, (r T - r 0 > 

(x T - x o ) 2 ♦ (y T - y o ) 2 

p m ^ t - x 0 ) 

( X T " X o )2 * (y T " y o )2 



take the 



(5.26) 



5.4.6 Weight Matrix . 

Since the direction observations to the ship from each occupied 
shore position are independent of one another, the weight array P is 
a diagonal matrix. Assuming z.iat all shore coordinates are free of 
error, for each observed direction to the ship the variance is 
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expressed as follows, as illustrated by the observation from station T: 

1 (TS ♦ 0 T ) e "Vs + "Oj. 



m. 



(S.27) 



The variance of the orientation, m^, with n known network point 
entation bearings, is given by: 



2 1 

m„ « 






2 12 

* * * * m T(n- 1) 



ori- 



CS. 28) 



(n-1) 2 

In this problem, two known points are observed from each occupied sta- 
tion. Hence the orientation variance is: 

(5.29) 



2 2 

m o T = \l 



But the variance of an observation to the ship and to known network 
points is the same for the instrument at any particular station. Thus, 

(5.30) 



2 2 2 
m TS * "Vl " "*T 



With (5*29) and (5.30) , the observation variance given in (5.27) be- 
comes : 



2 2 
m (TS ♦ 0 T ) * 2m T 



(5.31) 



Therefore, for the four occupied stations, each with two orienta- 
tion bearings, the weight matrix takes the following general form: 



X 

2 



1 

m^ 

0 

0 

0 



0 

™U 

0 

0 



0 

0 



0 

0 

0 

1_ 

"w 



(5.32) 



5.4.7 Normals Matrix. 



The normal equations coefficients array N is formed by combining 
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array A from Article 5.4.5 and P from Article 5.4.6 according to the 
following formula as explained in the previous section. 

N = A T PA (5.33) 



5.4.8 Constant Vector . 

The constant vector U, a column array, is formed as explained in 
Section 5.3 by the expression: 

U = A T PL (5.34) 

where A and P are developed in Articles 5.4.5 and 5.4.6. L, the alter- 
ations to the observed values of the observable quantities, is computed 
using (5.4). Its component vectors, L q and L^, are described in Arti- 
cles 5.4.3 and 5.4.4. 



5.4.9 Corrections to Unknown Parameters . 

The alterations Ax and Ay to the approximations to the ship coor- 
dinates are computed by equations (5.18) wherein N and U are as de- 
scribed in the previous two articles. 



5.4.10 Adjusted Values of Unknown Parameters. 

The adjusted x and y coordinates of the ship position are calcu- 
lated using (5.1) with quantities developed in Articles 5.4.1 and 
5.4.9, In matrix form, this computation is expressed as follows: 



X 




X ♦ AX 


a 




0 


T - 

P 

1 




y o ♦ Ay 



(5.35) 



5.4.11 Adjusted Values of Observations . 

The adjusted directions L to the ship position are computed with 
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(5.5). The matrix array for this computation takes the following form: 




(5.36) 



5.5 ERROR ANALYSIS. 

Using a procedure analagous to that described in Section 4.5, 
expressions for estimates of the variance-covariance matrices for the 
adjusted values of the unknown parameters and observed quantities are 
computed. These arrays indicate the accuracies of the adjusted direc- 
tion observations and ship coordinates. The matrix notations employed 
in the previous section are continued in this development. 



5.5.1 Variance of Unit Weight . 

The variance of unit weight is expressed as: 



T 

V PV 



(5.37) 



o degrees of freedom 
The weight matrix P is explained in Article 5.4.6. The vector of 
residuals of observations is computed using (5.9). The number of de- 
grees of freedom is simply the number of direction observations less 
the number of unknowns, the x and y coordinates of the ship. For this 
problem, then: 



degrees of freedom *4-2*2 



(5.38) 



5.5.2 Variance-Covariance Matrix of Unknown Parameters . 

The variance-covariance array of the ship coordinates is given by: 
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(5.39) 



_ 2 n 
E v s m Qv 
X o X 

where Q^, the weight coefficient matrix of the unknown parameters , is 
the inverse of the normals matrix N described in Article 5.4.7. 

5.5.3 Variance -Covariance Matrix of Observations . 

The weight coefficient matrix for the adjusted directions to the 
ship is expressed as: 

Q. = AN _1 A T (5.40) 

a 

where A and N are as described in Articles 5.4.5 and 5.4.7 respect- 
ively. Multiplying (5,40) by the variance of unit weight and simpli- 
fying the right side of the resultant expression by substituting (5.39), 
the estimate of the variance -covariance matrix of the adjusted obser- 
vations becomes: _ 

Z l a * ^A (5.41) 

5.6 NUMERICAL EXAMPLE . 

The solution for the determination of a ship position using four 
occupied shore stations is presented. Owing to the unavailability of 
actual survey data adaptable to this adjustment procedure, fictitious 
information has been generated for the problem. 

Coordinates of all shore positions were arbitrarily chosen on a 
plotting sheet. Observed directions were computed to the network sta- 
tions and the ship so as to permit the determination of realistic ori- 
enting angles. 

5.6.1 Problem Input Information . 

(a) Coordinates of occupied shore stations. 
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Data units 



meters 



(b) 



Station 


X 


Y 


T 


5000.0 


15000.0 


U 


14000.0 


17000.0 


V 


22500.0 


20000.0 


W 


30000.0 


23500.0 


Coordinates 


of other network 


stations used 


Data units: 


meters , 




Station 


X 


Y 


1 


2500.0 


4000.0 


2 


11000.0 


6000.0 


3 


21500.0 


9500.0 


4 


28000.0 


13000.0 


5 


32500.0 


16500.0 



(c) Directions observed from occupied shore stations. 

Source: generated using plotting sheet position coordinates. 



Observing 


To ship 


To orienting 


stations 




Station 


target 


First 


Second 


T 


34° 07' 50'.'5 


195° 58’ 20'.'2 


149° 28' 


40’.’4 


U 


15° 38' 40'.'3 


205° 22' 20'.'l 


145° 07' 


..M, 

00.3 


V 


346° 00 ' 0S'.'2 


200° 4l' 30'.'7 


157° OS' 


35'.'9 


W 


321° 54' 40’.'7 


210° 52' 0573 


180° 25' 
• • 


S0'.'2 



(d) Standard error of theodolite observation: 
Source: Manufacturer's handbook (Wild T2). 



0.5 



5.6.2 Auxiliary Information Computed by Program . 

(a) Approximate ship position coordinates (in meters). 

X Y 

15499.6 32500.1 
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(b) Orienting angles 



5.6. 

(a) 

(b) 



Station Orienting angle 
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3 Adjusted Values and Associated Standard Errors . 
Position of ship (in meters). 

X m Y m 

15499.84 0.12 32500.03 0.18 

Directions to ship position. 



Observing 


To ship 


m 


Station 


target 
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30° 57' 48'.'o 


l’.’6 


U 


5° 31 ' 36'.'9 


l’.'7 


V 


330° 45' 02V3 


l'.'6 


W 


301° 49' 38'.'0 


• 1 

1.8 
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6. SUMMARY AND CONCLUSIONS 



6.1 SUMMARY . 

Adjustment procedures have been applied to the two basic coastal 
hydrographic surveying problems; resection from the survey vessel and 
intersection from the shore. The adjustment procedures for each prob- 
lem have been developed using the specific observational information 
and characteristics of the particular survey situation. In the first 
instance a generalized least squares adjustment computation method was 
employed and in the second the method of variation of parameters with- 
out constraints was used. 

Background information concerning various aspects of coastal hy- 
drography was presented as a basis for the types of observational pro- 
cedures utilized. The various visual and electronic positioning sys- 
tems and associated accuracy determination information were described 
to justify the particular data input and equation developments used in 
the adjustment problems. 

Computer programs were written as the means for solving the 
adjustments. Information concerning the state of automation in coastal 
hydrography data acquisition and processing was presented to indicate 
how such problems could be solved by high-speed computers to provide 
the desired positional results to the hydrographer on a real-time 
basis. Modifications which would be necessary in order to utilize an 
existing automated survey system were outlined to indicate the feasi- 
bility of implementing these adjustment problems. 
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Certain other positioning systems; photogramme trie, satellite and 
laser; were described with their principal capabilities and limitations 
as applied to coastal hydrography. With the improvements being made in 
these and other positioning systems, it is envisioned that such me- 
thods may in the future replace or supplement the present day electro- 
nic and visual systems as the input data sources for problems such as 
are dealt with in this investigation. 

Finally, numerical examples were worked using computer programs 
to indicate how the adjustment procedures could be applied to actual 
survey situations* Because real observational data was unavailable, in- 
formation was generated for the purposes of these problems in as real- 
istic a manner as possible. Production of this data included the use 
of plotting sheets and instruments utilized in practice in conventional 
coastal surveys. Arbitrary instrument and system accuracy estimates 
were used in the numerical examples, precluding their value in giving 
any more than a general indication of the accuracies which could be 
expected with such adjustment computation procedures. 



6.2 RECOMMENDATIONS FOR FURTHER STUDY. 

(a) Observational data from a single sounding line was utilized 
for the Chapter 4 adjustment computation. In practice, when sounding 
information is required, a series of such lines are usually run to ob- 
tain a satisfactory coverage of the area with depth determinations. 
Observations to shore points and electronic positioning data obtained 
during these additional runs along the same coastal region could be 
sequentially combined with the first group of data to strengthen the 
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results obtained for the shore positions. 

(b) For the Chapter 4 adjustment procedure, it was assumed that 
no positional information was known concerning the prominent points 
along the shoreline used for sextant angles and compass bearings. Very 
likely, information would be available concerning at least a limited 
number of the terrestrial points, either as the result of earlier sur- 
veys or from aerial photography. The inclusion of such additional 
position information, provided it were at least as accurate as the data 
observed during the survey, would serve to improve the overall posi- 
tional results determined through the adjustment. 

(c) In the manner outlined in Chapter 4, the ship problem could 
be used to ascertain errors in electronic system lines of position 
caused by variations in electrical ground conductivity along coastal 
areas. For such determinations, all observations would be taken to 
geodetically determined points ashore. With the known shore positions, 
the need for compass directions, the weakest part of the visual obser- 
vation procedure, would be eliminated. Such a procedure for deter- 
mining fixed electronic system errors is given in (21). The method 
described therein utilizes only individual sextant angle positions, 
however, with no overall adjustment of data. This problem applied to 
their method would certainly appear to be a means of improving the 
positional information used to determine electronic system position 
line irregularities caused by ground conductivity variations. 

(d) The use of photography of the coast taken from the ship 
would appear to be a means of improving solutions for the determination 
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of ship positions such as dealt with in Chapter 4. As mentioned in 
Chapter 2, terrestrial photogrammetry has been attempted on occasions 
with little success. However, with sufficient emphasis and attention 
devoted to developing a workable procedure, the use of photogrammetric 
information could enhance coastal survey positioning. Real-time film 
developing techniques and shipboard computers for data processing 
should allow photogrammetric methods to be integrated with the conven- 
tional procedures now employed in coastal survey operations. 



6.3 CONCLUSIONS . 

For coastal hydrographic survey situations wherein shore control 
is unavailable, the Chapter 4 adjustment procedure affords the hydrog- 
rapher a means of obtaining both terrestrial and offshore positioning 
information. This adjustment method is primarily intended for use in 
a rapid, self-contained ship positioning system which could provide 
satisfactory positioning data along an inaccessible coast. 

The actual position accuracies which would be obtained with this 
procedure must await evaluation of the method in a real survey situa- 
tion. The positioning improvements may not be great enough to allow 
the results to compare with the accuracies requirements of conven- 
tional surveys conducted with the normal shore control. However, the 
accuracies will quite likely be better than the results which are now 
obtained using either visual or electronic methods separately in such 
situations. 

When shore positions are known to geodetic accuracies, both 
adjustment procedures developed in this investigation can form the 
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bases for ship position determinations considerably better than are now 
obtained. While the present methods provide sufficient accuracies for 
nautical charting purposes, they are not satisfactory for certain pre- 
cise offshore positioning requirements such as for engineering projects 
or emplacement of scientific instruments on the ocean floor. 

Both adjustment procedures afford a means for the initial cali- 
bration and later monitoring of electronic positioning systems. With 
sufficient shore control, the offshore locations determined through 
these adjustments can also be used in analyzing irregularities in elec- 
tronic position lines caused by ground conductivity variations as the 
radio waves propagate over adjacent terrestrial regions and then the 
coastal ocean surface. 

As described earlier in this investigation, random and systematic 
errors are present in both the visual and electronic methods of posi- 
tioning which render the Chapter 4 self-contained shipboard procedure 
unsuitable for more than approximate positioning sufficient to satisfy 
operational requirements. However, continued improvements are being 
made in electronic systems, particularly as concerns the determination 
of and correction for radio wave propagation irregularities which are 
now one of the greatest causes of positioning uncertainty. Also, other 
systems such as lasers and satellite navigation methods show promise 
of one day being capable of accuracies and characteristics which will 
render them useful in coastal survey operations. With improvements in 
existing methods and the anticipated new systems, there appears the 
strong possibility that these or similar adjustment procedures could 
in the future be utilized to provide the desired position information 
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with accuracies which meet the specifications required of the normal 
coastal hydrographic survey. 

Adaptation of both adjustment procedures to «<»rvey situations 
wherein the coastal positions are accurately known will aid in pro- 
viding offshore locations with considerably better accuracy than can 
now be accomplished with existing survey methods. While such increased 
accuracies are not required for the usual application of hydrography, 
the construction of nautical charts, they are required for certain 
other applications. With the greater attention being paid to all as- 
pects of the use of the ocean and the coastal sea floor, numerous other 
precise positioning requirements not yet dreamed of may well arise 
which will necessitate the greater accuracies which these determina- 
tions can provide . 
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